{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Elevation profiles extracted from SRTM over Mt Baker compared with ICESat-2"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 53,
   "metadata": {},
   "outputs": [],
   "source": [
    "import numpy as np\n",
    "import geopandas as gpd\n",
    "import rasterio\n",
    "import matplotlib.pyplot as plt\n",
    "import glob\n",
    "from topolib import gda_lib\n",
    "import topolib\n",
    "import gdal\n",
    "from rasterio.plot import show"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Read reference DEM"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 25,
   "metadata": {},
   "outputs": [],
   "source": [
    "dem_fn = '/home/jovyan/data/srtm_elevation/SRTM3/cache/srtm_wa_subset.vrt'\n",
    "dem = rasterio.open(dem_fn)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Read ATL06 data"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 40,
   "metadata": {},
   "outputs": [],
   "source": [
    "data_dir = '/home/jovyan/data/nsidc/**/'\n",
    "ATL06_list = sorted(glob.glob(data_dir + \"*.h5\"))"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Sample DEM at ICESat-2 ATL06"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 41,
   "metadata": {},
   "outputs": [],
   "source": [
    "ATL06_fn = ATL06_list[0]\n",
    "dataset_dict={'land_ice_segments':['h_li',\n",
    "                                   'delta_time',\n",
    "                                   'longitude',\n",
    "                                   'latitude'],\n",
    "              'land_ice_segments/ground_track':['x_atc']}\n",
    "\n",
    "ATL06_gdf = gda_lib.ATL06_2_gdf(ATL06_fn,dataset_dict)\n",
    "ATL06_gdf = ATL06_gdf.to_crs(dem.crs)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 42,
   "metadata": {},
   "outputs": [],
   "source": [
    "poly = gda_lib.dem2polygon(dem_fn)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 43,
   "metadata": {},
   "outputs": [],
   "source": [
    "ATL06_gdf_sub = gpd.sjoin(ATL06_gdf, poly, op='intersects', how='inner')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 44,
   "metadata": {},
   "outputs": [],
   "source": [
    "points_xy = list(zip(ATL06_gdf_sub.geometry.x, ATL06_gdf_sub.geometry.y))\n",
    "\n",
    "refdem_sample = []\n",
    "for val in dem.sample(points_xy):\n",
    "    refdem_sample.append(val[0])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 57,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAfAAAAHwCAYAAABZrD3mAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjAsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+17YcXAAAgAElEQVR4nOy9e5Akx33n983qxzz2vTv7fi8Wu0sACyyAFSkRkkCC5J0IkpBph3Wkgz757BB5lGWLjJAZpi+sOIbDto6k4siw76TjQbo/zDuKxzApnSjq9DAPlCCJohaLJ7HAAvt+zs7s7MzuzqNflf6janq6a/KbM9nb0zM98/1ETERnVlVmVnVN/Tqzvr/fz1hrIYQQQojuIlrsAQghhBAiHBlwIYQQoguRARdCCCG6EBlwIYQQoguRARdCCCG6EBlwIYQQogvJL/YAQthYLNid/T2z6k0U7gpnDDmG1NP9FxNr3PWk2oeJYlJPzttzOdgmOiy6gX1HwFW7Cbewpl61BuPYY4bq39/FeAvuYFXT9r3RYEDfbKyeE2dthdbnyIbI83vbkG3smCjnrLaG9e3eP+mbtRXWtzH8cWRIH+y8TeAXa+G+/wGAutraKtmftMXqPf9MJq6RQ0g92Z+2E/PzBjtvdgytD3+G0K+DPPNsTL5v9owEYNk20hbbn7bTQt+u+isTJYyUKs4DusqA7+zvwbefPDarvthbcu6fI0YJAKKce1uuWHHX593/rK38eAiF3ZxxjTwIPefNKK6adNbnesvuMZG+AT5ekyMPF3YNSX2Uq+GleD/+Ye1/RhXJOPKo4f/O/QYe7XkTAPBPSx/DN2tPIbGSFv+g8Of4Qt+/cYyJXCvSt/fa5skxPeTHUQ/5B189+0cqAKCvn3Zte/vc9f2r3PW97vq42Ouu713N++5ZS9py94Hiemd11LuZ9pEvuI/JRe5rFXl+DLioxe5nSLJtyllfLY866231LmnI/T+Gqrt9AIhKY856U55w7z/pHlM0cZvUk7ECMGVyTSbdfWOCnN+k+5lq+WnDltw/zGyp4KyPSX2tVKR90GPKrC13fZXsH1f5PVirup+fVccxz/x/r9F2tIQuupJj0Rk8iZfSkkEVOfxB/O769ndE59NPiVF9oF4WQojlgQy4WDYM25mZ4Kl4X/opmeX+RfWRzg9ICCEWEBlw0bVsMre95Uaeqz2GF6v3L/SQhBCiY8iAi67lAXMh/WQzZeDn888jh1q6zSCGwY9q7+j4GIUQYqHoKhGbtcb5kt+U3OIhW3ALzwAgD/c2Js3K9biFGK2o06lqMZC8cQvMmFgMAPL9buVIRM4vt8q9v63wW8cSgQZVitKGmMo+uean4j3TFQAsTpk99fN4HKfwj/BdPFv6eQAWFhE2Fkbpec6CCdJ83zc7hl2qXrKhSIQ3PhU620aUwaZM7gPeA8VShbO73sbu/72YqLoBoEyEYVHeLa6LIrcYz0Tua84EaQBgS0PuDUx8Rs7DkPM2VS6gY2I1UyEitqlx0g7po8qvOcru5wutr5Lvm+xuK/xuY2K12jgRWZLnEROeAVxkxo6pkf1pO+Q5CACVCuvDcYznkaMZuFg2NL4DB4C7dloFnfwQ+HF1f4dHJIQQC4cMuOhafj7/PPL1ZXLgB/GjTe+5h+N1Tftny0II0c3IgIuu5dHc23gyejEtJa5kv1/62fr2gah5WXQgcvvUCiFENyIDLpYVjbPsB/Ln00/JDH2NIQEohBCiC5EBF13NQMZ1rHHWPRqvgUGM6Xfgvzv5EbxYOdTJ4QkhxILRVSr0dsJCqUZEwc3qaSjOFggNy2qIyj7XQ2SfnmPye4katUzOr8KVszTkM+mCqVFt2X17Rn0zSvKj46fxzbGnMD3Lfqj/bUSrkvP/yfzLiKZ+AbXUlayGHH6//DN4rP8NOvaZTlh4V34IVZuz2OYMpvL1xK02FeKRUHAr2m3RHYLUFN3fqynyuJdxL1H2MxU6UV2zegCIi24VekwuetyzgTTkHmt09wrtm6nEaXzxdsIU7VV2j5DnVJV40UyR8KcAQO4plMk/8qT7/rTjbjU2C2UKcC8XFj6aqdB94Uypcp2pzSvuvkPCos51TM1xftYT118zcNHV/LgyrSxPleaVA/Vtj/acxlM9f7cIoxJCiIVHBlx0NcO19d7yk30n00/JjPrBwrlODEsIIRYcGXDR1QzkeAAOoHFGPj1Dly+4EGJ5IAMuupqP9v8gjaqXzLD/09TjeLE0I1Sba4YuhBDdigy46GoeLZ7Ge3teSEsGVeTxnYn31LdnfcGFEGK5IAMulh2Ns+yPrnqueYZeehwvluVKJoTofrrKjcwYi0JhtjtELu92nSh63KlyxbDkJMx9odBPXF88iTtCk3qYPHEfIu3YmP8uY+5z9iZxP2HRRz13Dj07egxJuFEgrlGZfAYDd0aBhq9hIDda/2n6aN9pvLfvBfzZ5LswM0N/Eo8WTwMAbI2MliRSMTF39bPE9czUyDGTxP0qR9x0In4/U4puVx3qFtbjThYB4nYGADmSKMOU3S5KtneVsz4mCTqSttzbbLHf3RZxv4omRpz1ufGbtO+IJH5hLluISDIfVp/n7lSsLeoWxr6LCXeSE0x6ghsxd7ESSUYzSdxBibuWLyFSTFxIbZX0USP1noQiNeZ6RtzF4nj+rl8AEJMxAYC17m2x69mtZCZiOfNgz9n0U6o0L57lO0PvwYUQywMZcNH1/LjcrDT/weRjTdtnxUSfQ7kuhBDdgAy46H4yS0zfnzrepESfmZHLF1wIsXyQARddz0fXPodcPa2oQQyDv516qL496wv+g6lHOz5GIYRoNzLgout5rO80/rt1f5CWLCwibIgakpxkZ+il41KiCyG6nq5SoUe5GP2rZ6smDUk0wtTpAJAnKnQaZJ4E3o8ra531TM0OcHViRBKN0MQrZH92PQBPkgCioI7GiBLVd+ewnBtMNU8U3wzrENTerqxOPxkANnkvno7jo/3P4d+Pvx815DA9Q//hxMN4xJ4FmFKUfH9MaQ6gSQk/Hwy5P02BXMA879uwn+Is6UyZKOCr5J5iiS0AgCifI5bMhCRlYcpqALB5omgv3XH3TdTb+dEb7g48fTNldzARUUp7FP5Mhc7GG024k75QtfmU53utEI8VJlyvMq+YwGQ+HYI9d5hCvBPn4RqTkpmIZc9wbV2mPKM0f7TnNP5Rz3fTUjJDX4/mNKRCCNFtyICLZcFAbqy5nFGe38W0v3Dya/b1mmKiCyG6GxlwsSx4MN+sNJ+wzcuSw/E6b1kIIboNGXCxLLhl12JahQ4AfzD5Hnzz7vvr2wei7Ay9uSyEEN3GvA24MSZnjHnRGPPdtHzMGPNDY8xLxpgTxph3Oo45nG6f/rttjPlMuu2bDfXnjTEvte+0xErjXcUfpzfzjBH/j5M/Vd/+QG7a9zuZoa+GJ4SkEEJ0ASEq9F8FcArAtOz6iwC+YK39Y2PM02n5PY0HWGvfBHAMSH4AALgC4Dvptn8wvZ8x5jcBzDklMsYi74hvniNqbF8s2tKEO+azJXFtQ6l44vzmifo4yoep0JmSPvbFGCbbcj2kLaZ09yj8DVFqUxWnJ3b7fHkEF/Hf5r+LZ6sfwbSRfn98EtWRJO72yNQmmFTABlj8m9JH8FT1NRzLv03GRMbqUaGbiKirWfx5EnPZVIk63XPNwZTrpNpQVwGiSiZjAgAUyT3C9icqdJ8SnMViN3lyDItTzmKIszH5YGp61jfbf8qtsAdAlev0GpbJ98fU5kRpDgCWCdRbyMEQClWCk+cz69vn4RKqKueq9XDvmpj0HQd65MzrihtjdgH4EIBnG6otZoz5OgBX52jmfQDOWGsvZNo2AH4BwDfmMxYhGO8rnKwHdMmhhkPR5fq2d5o3EMFieoZeQ4Q/iN+9WEMVQoh7Zr4/mb4C4HNo/j3/GQBfMsZcAvBlAJ+fo42PwW2kfwbAoLX2rXmORQgnv1/56dTXO0INOfxB9afr245FZ/AkXmza3+PRLYQQS545Dbgx5sMAblhrX8hs+jSAz1prdwP4LIDf8bRRBPAMgG85Nn8cntm3MeaT6Tv2EzenPMt4YsUzbDNKc9scZOdno1fTT4npfsA0LQYJIURXMZ934E8AeCZ9z90LYK0x5usAPoLkvTiQGOZnyfEA8EEAJ621g42Vxpg8gP8cwOPsQGvt1wB8DQCODfRp0iQoA8Yvozhl96afkmhtM2UhhOg+5pyBW2s/b63dZa3dh2QZ/PvW2k8geef9ZLrbUwB8S+Bslv1+AG9Yay87tgkRxH9WeB75elIT4Afxo3ixdrC+PTsjz5aFEKKbuJdY6L8E4KvpLHoKwCcBwBizA8Cz1tqn03I/gA8A+JSjDfZe3ImFW1leLfXP3hk8rjkAVKskJjhRrjNFIWsn8qiVy4bENicqZqY2z5eKzvpcgat5c5NuNW+O9MHqo3wLqt0F5iFcx8+aV/F9cwwwBlWbw3cm34Oj9hoAYJMZTybfyQQcm+w4apPua8hgCnsAVKEesdj0bP+YxJ/3qmaJEjxmHgHueqrGrnrOm7TFYH14Y6GTbTbvzlGAwHrWPgCugCfnQa8GUad7VegkNj1VpzPyLKa677Vk4PdKnl80BryvLfKM9P7/OdvxeY2Q/wHSNxDmoeRVwLegXHcRZMCttc8BeC79/DwcS9/W2qsAnm4oTwDYRNr7b0L6FyKUm5iZZR+xF1PjnfzjrrITHl8nIYRY2igSm1hWDGSSlDSWx8xqGFjAJFb7/4k+gJdwAEII0Y3IgItlxRF7MfmQzrLrZQDH7WnkECfbTJJW9ASUF1wI0Z3IgItlRTLLjgFjECHGmFld33YM5/C/xN9AHjVEtoYiajiO04s4WiGEaJ17EbEJseQ4bk+jx9RQsUABNRy3zQb6v8Rf4f74Kk6YQzhuT+OYOUtaEkKIpU1XGXAbR6g6lNflsltJXK0QJSp4rPJaYExdFrPXB1M55olaOUfOI09ip+cL/LxzpA8Wb52p8ntWceUsVaMSqlOBSnAyVgB4oHcQv5X7LZzIH8Tx6tt4oDaIMvphGs77QVxFNSrib/MPILZ5PBKfD+ibq3aZMp/FPWZtWaYw9imPy+RfmcRPZ7HQLdEGm6JH/cuUzGWPsjuQUK0hU6fbojsHAlWtAzyuOlOhhyrEV/PUtjbv/t+IJm47683EuLuedeBRoZvI/XxhdwKNRx4YcxwIf66GPnMAv0K9Hfv7CFWbM7rKgAsxHx6uJRHWTuQPNpWneTnah3/c/9+jghwK+Pv4Wun/CjLiQgixFJABF8uOV3J78anVv5wa6Br+1d1/iUcws1T+Qv4gysjDmggVC5yI7pcBF0J0HRKxiWXHd4vH6wa6jDy+WzzetH1dPA4LA1iLGBHW27uLNFIhhGgdGXCx4ngjtyv5kPqDn4p2L+JohBCiNWTAxbLjSDUNrZ/6gq+OmwV3N80ab1kIIboBGXCx7BiNVgMNEde+3vtevBztq2/fZO807Z8tCyFEN9BVIrZqNY/BwS2z6iu1MJcw37aYyPtzgW4KkcflgLUVEfcytj91I/O4OjFXiDxzNyJjqnhc9HJkXAx2rdhYWYKOaY7VLiDXa1FLI67VbIQ/zL0LD01dAQAcMteAAuoz9EPla7PcE5lbnfG6f7jPm96FJJkCc8eBJzkPS8xg2TEskQNp33qSSNAEE8xFiblZ+e6bqvsY+m1EgYknfK5fxJUrZn2QtuJin7tv0j4AGObCRgidkRnre65NuKtJsh1TIW51MUsgw78jS77vmNXX3G0xN1gg3JUrdH9vIpU2uaRpBi6WHUcrF/HE1Knmyob/lzfzO5MP6Qz9rwtHOjQyIYRoHzLgYlny7vKbyYd0ln24dmVmY+bH718UHsAr+b0dGpkQQrQHGXCxLMnOsutlAB8qv9CU1MQiwgv5+xZjmEII0TIy4GJZcjNa3VxuSGrycPUC/qvJ55KCtbAwWGfdISiFEGKpIgMuliWbYn9wlvEoFRRNz9CnfcOFEKJL6CoVemwNStXZisYqUSD6VOhMURgTXWs1dvcR0dD+HKZAZGpzVm8c1wIAckQ5DngSphDVdS7nVpwWalwxHBGFOlO00+tBlPFMfZ+0ldzSPzf+Mv5D/3HUbPK9PV98AC/hAI5Wkvzgw3Zt03HDZg1Vsjbi/cXrUf8HwZI/RFwFyxS9hlxD2keuBXUsu92YtwBTp1fDvBeSPsg9RepRLrn3z/NHIUuMYojanO0fseQnxfD7hqnTTeg1zHlMQK9bHc/vQqJOZ2MiCaUAwJL/RVpP7n+fcpwp1GPyrK+xvtuUmKQVNAMXy5Kj5Ut4YjIVshmDGnL4Xu9j9e3ZGfqmmsKpCiG6CxlwsWJofC9+uJqq0lOV+irL06MKIcRSRAZcLFs2ZmfZDeWxaBUao7X92/4n8Up+TyeHJ4QQ94QMuFi2HCpfSz5M+4JXZ3zBHyufRQ627koWI8L3eh9fjGEKIURLyICLZcvp4vbkQzrL/puew/Vtc0VrE0KIpU5XqdCtNSg5FM5MIV5jStsOYDxds9jfTM3I1PRMnV4zHvU90ZBa61a1xuQa+pSXNNY7U9mzuOM193WqGa4Wb4xhnh3jXxaP4JX83roS/adKp/EXvQ82zdCnj2Hn51e1knGR7zsi/35M1B154pGzPqjanHlokGtuqrxvtsXQLUR1PVWmfaBIrkpEFM4kPwJVRHtU6CZyK9cZlsVIz7uFkra3n7fFFPBT7jjlrB5lcm1LU7RvTJBtJOY5wxSJOr3E1fdRwX2MJcpx6lXk8SxhivZqoKI9VM3uaysUzcDFsuXpiZcQNUZcMwYni/vr28eifhjEgDGIEGM0WrWIoxVCiDBkwMWy5aHyJfzarT9EztZgbA1FW8Nj5XP17Y+Vz6GIGnK2hgJqeLx8dhFHK4QQYXTVEroQofz8+AsAgOf6H8B7p16rL58DyXvwf3HrX+N7vY/p9bcQouuQARfLmteKu/HV9U+jYnJ4ubgP91UHm4w4APxh33FUkcN3+47jt259DUerF0lrQgixdNASuljWnOzZj7LJIzY5VEy+6R04AHyv9zFUkIc1ESrIN0VrE0KIpUxXzcAtDGp29m+OmKm3Hfu2G6a09cZIZ8JgsoGJj0PV6QC/VqyeXUGvGpttI+NlsdAjEvvb17fNqK7XVScS5b21iGFwF71NyuxZWcvSsiVjtUztDU88ZnIMVZsTL4LY8+8atcvjgqrWPfG6iXIdVXKGRNHO1MpJ/0yFTu5QFgM+Impsjwqd9kEwzAuEtGNIfHYAsMWesGOmiHKcqc19yv8Sid1OmqIK8dA4/QByfeT8fJ4Yrr49/xe1slvhT9sK9ExhHjy+Y0LRDFwsa8Zy/WiMuPbv1vwMXi3urm/PRmvLloUQYqkiAy6WNY9OnU9WQxoirv1x36P17UcyMdGPNERrE0KIpYwMuFjWPFS+hCcm3myqa1yEezO/M/mQztD/ungYQgjRDciAi2XPT02+lXyYjrhWuVrfln2j9pc978CrSmoihOgCZMDFsicbE/10YUd924emTiLXGK0NEV4oHliMYQohRBBdpUKHdSuvWcxzqoZuARa/vBVCFd8szjVTb7dCqMKy5osxTNT/NnKrTlmM9FZwjfdmLqM0N2vrcYofnLyCj+f/El9f8yRgLSwM1lYmW4qFHhormXkLsBjN8HhVsCtoyP+GIX0YMiZD4tUDXE1sWYz7IlEfM9U6AMNioedJfZmp0EkHeY8CnnhD8P2ZCp2044lHbnzqeBdTTJ1O1OYVj1cFU5tXyP93ya3qjkk9U6cDQFxxn3eN9UHvZ35++R6PAt8B80ypISx2OuDxBnLET/e1oxm4WPZsrI17t49HvckHxwxdCCGWKjLgYtnzc+MvIYda/R34X/UfwmsNrmTMF1wIIZYyMuBi2fNQ6TKemDidFIxBDTn8cf+x+nb5ggshuhEZcLEiGWmYZddV6ekMfXXsyZEshBBLBBlwsSLYkJlVN5bHolWYFa1NrmRCiCWODLhYEWRn2Y2+4I+VzyKXidampCZCiKVOV7mRWQBVh5S/ne5iDOr6RbsOd/FibmE02UcH3MtYwhS/W53btScwJ0RLrlzsezqdT5XlJkls8sOeQ3jmzklYE+Gh0hU8MfkG/qLvgaY+XIkQLEtUAV+iE7ebCUvlEJHEHSbHkz+Y2O1ew9xoTI25WZF7zds3cVVjySpY0gufqxoZL3VvY+5D5OszPjcyBrsVTBufR6wP9uAhiWJsiVyPssctkrmLMRevSXfildpk0VnPXMUA7krJjqmUWB/cVa1WdbdVJfUxS3BE3DurHlfbWo1cQ/K/wdAMXKwIso+v5/ualeg/OZWJ1qaY6EKIJY4MuFgRfHD8ZUSNEddMhBd799W3Z6O11WOkCyHEEkUGXKwIHipfwsfH/ioppBHX1tUm6ttH5AsuhOgyZMDFimFWxLXpWTfk+y2E6D5kwMWKYSQTE72x/MGJTLS23iN4tSBXMiHE0qW7VOjWoOIM9u7ev51iUIZlnXt+GuWISpuqrhGWrCXXggqd9cHE9L70I1ycGzYupu6Ma7wdY9zKz8hYbMjERN9QG68nGnlw8gqemDiNv+h/RxKtzebwRz2P48HJq5mWwtXKxhJVPtmfJibxqO+pBwNLFFN1985U6zbHb2iWGIUp15k63RRaUYKTa+JJYuHCuzdTrof+nwWOyQ/7xyTPEKb8L/uU4OT/jyi+qdq87PaQqHlU6ExtXiVtxUQ5XiZjBYAaSabCFOJVsj9TmzM1O8Cf9RVH30pmIgSAQ+VmX/B6mZB9Ly6EEEsJGXCxYpiVF7zhHTgAbKgpJroQonuQARcriOxSVHP5g+MvI2drgI2RQw0fnHipc0MTQohAZMDFiuHnxl9CwdZgbIwIsXMJPUrfLVoYnMlv7fQQhRBi3siAixXDQ6XL+MzIH8FYixgGX9n4oaZobC/27kPV5AATIUaE31z/YbzasF0IIZYS3aVCB1AjKkt6wALDBII+xXCuA+p4BlM0MuUzE/maFi4u65upONn+EQ9APydvFHYiNhFgDCow+N6qY3jHVBI29eHJC4jWW9QakpqcLBzAg5ONYVV9ylK3ujoiX7glsZUNUf/6VM805rknhnlIO5FXhU68KsgxVLVOrkeyza1Qt0RlT9X07Pnh6dszqKDdadz2VmAeK4Hnx5TmAGCJ6popxFk9U5tXp9yx0wGPQpyo0Jniu+pRulcq7rbKpJ71wfJF+NTjNeJh42qLOToBmoGLFcZIbhUtP1S6jF8Y/euk4IjWJoQQSwkZcLGi2JjxBc8ynuPR2oQQYikhAy5WFH//ziup0jxZl/qbvkN4rWdXffvsGbp8wYUQSxMZcLGieHDqCn5yIk0dagxqJoc/Wf1IfXt2hp71DRdCiKWCDLhY8TTOuu8vXUs+pDP01fHUYgxJCCHmpKtU6IBb2RcHxgpvhYgogD1huYPbYjGzQ1XrPvUjVZsH7m9Y7HQAMVFlMhVnRBTDcezuPR+orAaar8mG8iTQsFK+oTKBWjq2MbMagAVMBFiLb6x9Au+++zYeLF0GgGT5nRAxZTC5SbhCnMUW5ypmGgs9UNEerN4GEJHvyZBjInKfx54+bGDsdhp3nMUKD/FuaTP0HAB6HkzxTWnlvJlSukRimxOFeGWKxE4nam/AE3ecxVUn8csrnutUKrvHVWFqc3I9mBeNz/6wZ7Srj9gzz9YMXKw4DpauJx/SWfb902UAj0xeSIK5TLuSmQh/subhxRimEEJ4kQEXK463e7clH1Kl+Vs92+rbHpy6gp+6ezpzRAcCCgghRCAy4GLFccvjCw4A75o4k3xwzNCFEGKpIAMuVhwbqn5f8PqMvD5Dly+4EGLpIQMuVhwfuP1aky/4D1fdjx/37qxv90VrE0KIpcK8pYzGmByAEwCuWGs/bIw5BuC3AfQCqAL4ZWvtjzLHHAbwzYaqAwB+3Vr7lXT7/wDgV9Lj/8ha+znfGCwMqg6VXo2pfz0qQF9caResjxxpx6dAZLFzc4aoc4nim8Vb90ZcZteKqZLZdfL89GPKVsNif5M+coGxpgF+3RvvhUN3BvHOO2fxN2vvT3zBkfiCH7mbLJVvqEw2Hxubemxm3z0VEYU6Ow9fvHznOXhcHgy5d1j8ba/y2YEvFjry7vOm58eU/ERpDgBR3q30pZDrQSEKY4Dfz6HKdUvUyt5j2P1MVNqUwJjxAD8/Go+cqM3LRO3NFOVAuKqcKcFZvHMAqMYk3jp5TlVIH6HtADy+uUuF3q5Y6L8K4FRD+YsAvmCtPQbg19NypmP7prX2WLrP4wAmAHwHAIwx7wXw8wAettY+CODLAWMRoq00vhefNUNffbBphi6EEEuBeRlwY8wuAB8C8GxDtQWwNv28DsDs5MrNvA/AGWvthbT8aQC/Ya0tAYC19sZ8By3EvZJ9D95YfmDyCn7y7ttJIY3W9qdrj3ZyeEIIMSfznYF/BcDn0Lw6+xkAXzLGXEIye/78HG18DMA3GsqHAPyMMeZvjTE/MMb8hOsgY8wnjTEnjDEn7lQVFUu0h4OT6e/FdJZ9cA6leVa5LoQQi82cBtwY82EAN6y1L2Q2fRrAZ621uwF8FsDveNooAngGwLcaqvMANgD4SQD/E4B/b4yZ9dLFWvs1a+1xa+3xNfneuYYrxLx4u29r8iG95eq+4SkbZsVE9yvXhRCi08xnBv4EgGeMMecB/B6Ap4wxXwfwiwC+ne7zLQDv9LTxQQAnrbWDDXWXAXzbJvwIyex+IHD8QrSI9Rb/3u1X0/fgMXK2hr93+9XODU0IIebBnAbcWvt5a+0ua+0+JMvg37fWfgLJO+8n092eAvCWp5mPo3n5HAB+Pz0OxphDAIoAhoNGL0SLvH/0dRRsDcbGiBA7l9AjAAbytRRCLE3uJZnJLwH4qjEmD2AKwCcBwBizA8Cz1tqn03I/gA8A+FTm+N8F8LvGmNcAlAH8orU+wTwAC1Qd0nzmMOJz+WHuVKHUiAuUa5z1Y0gSi1DXNpLzAjmPIxl1Swl1q/OcH0vWEuoWFpOkF624B9pCpal88PYNfPLSf8Jv734KsTH4l1s/gN3jt/COiSQb2XlTsI0AACAASURBVMme/aiZCNZEqMHiZM9+HLpzA7nYk8yEJA6JPUlIQs4h50niYoy7b5Y4hLmd0SQgLUCTuzA3Q1/fniQybsLcrKzX5SfsPJi7WMzuD497GTsmlNDkNYDHjYy4ZpWmepz1VbI/S1gCeBKKMFcu5kZG3NEA/gyj9TSZCXE7I+5lvrZcFtD3vAsy4Nba5wA8l35+HolrWHafqwCebihPANjk2K8M4BMh/QvRTm4X+mCNgTURqrB4ZdXuugF/ePwSCraGCiwKNsbD45cWebRCCNFM16UTFaJdHL1zGfltNVRhkc8Y6Qcmr+L/PPdN/Pn6BxdxhEIIwZEBFyuWIxPX8Kmr38dfrT2EJ26frs++G/mzDUdRMTn82Yaj+Gfnfg9Hy5qJCyGWBjLgYsXyRv92/KsdT6Ficnh59R4AwAdHZtTmf77+QVRMDjAGFeTw5+sfxNEbMuBCiKWBBLZixfLqml2omFwiVDMR/sWO9+NUvzKPCSG6g66agccwKHnUou2AKf7C02pwQs/B5ogalCiMWVIUAMgR1WmFJJJgivLZIXfmPoaNKyLqdKoGjfj1Y2pbV1uHRwZhtgPWWsAYxIjwp2uP4r6xmyjkKzgwMQRsRF0aemBiCNaT9CIml53dU8EJdTyq3YglLaky9TG5d5inQL5K++bJPsiYiGqXnQPQQvKOQHxKX3Z+TD3OlONUze5RoXs9aVz7s++iTV43AE8QUi65k5ZwRblHIU6uIVObs4QiLAEJAMSBavMaubZMbc7a8bXlTGZCklkBmoGLFcyhu4P4idGzzZUN/ytn+rakdUnliTX7OzQyIYSYGxlwsaJ5bCzNrTM9yx4fovv+aO0BvN63oxPDEkKIOZEBFyuas/2bkw/pLLteBvC+W68ngVCml9iNwcv9exZjmEIIMQsZcLGiGc330/I7Jq7ho0MnkoK1AAzW1iY6ODohhODIgIsVzfo5DPJELg0PSbKWCSHEYtFVKnTAraastVFhGTOlaOD+rajWS0QxuTYTx3uaPFNve65HITAeOVOOMzW795jAtgqRO/41u+YAV8Aznhw6jT/f9ABqaezsE+v34fW+HThavghgdh7wkdwqqpB1ZMP1EpHzZnHKfQr4GgkVztpiyn8Wt92nYs6ROOWW1LO47XGNf3ehin2msmcxz33nFxyPnCrB29c3fU4Fnp/v3OgxTHXNVPlM7e3pO1RtzjxWWN+ARz1Ovj/aN7lOvnwRVaZ0d7TlyxCiGbhY0Ry+ex2Pj6ZCNmNQMzk8t+lwffuGaiYveFV5wYUQSwMZcCEyjBZm3oMfmLyRfEh/Bq+qlRZjSEIIMQsZcLHiWV+ZoOU7+T4Atv4O/P8d+Amc6pcrmRBi8ZEBFyueAxOp73dDxLVpjt69nPyT1F3JImUoE0IsCWTAxYon6wt+ct3e+rYjE9fwztG3M0eECaqEEGIh6CoVukWYutunxq6GxhhmSkOm1PQ840MV6mysPURp61PsMmUki4FdIG0Vc0T2DK50j4jhyxPlcy1yn3cxz/vOE+U6jUNtzaxRnVi/D68P7cSR8SS96GO3L+KH6++vz9D3jw85VbI8bry7nsVOZyrt1iDfN1Wb82sbSmjMc6aYTw4K/NEUhz3afAr/0DjiocpxH+1SlfvufwZTlbN6phwP3R8Ij23O1OZlTx8sVjn7nkLV5hWPCp327ajz3fmagYsVz3uG35wVce211bvq28/2N8dEP9u3dTGGKYQQTciAixXP4buD+Mi1F5NCGnFtTXWyvv1WJlpbtiyEEIuBDLgQACbzzRHX6rNuABuqCp8qhFh6yIALgWbfb6B5lv3ekVNJtLH0HfjfrduPN/q3d3R8QgiRRQZcCADrMr7gjbPuI+PXcHzsfFJIo7V9f8MDHRydEELMpqtU6Hljsb5YnlU/PNXr3L8V1SdXGgYqE70q9DDlepmoGZlCvBUKROWbI330ea5tL1Gos7bY9WCq9ZxH4s+OiVns7/Q89t8dBragPsvurZSa1bOZ8x3Jr5qlfLWh8boD9/fFeWf3OjuGqY/z7IlQ5eNibbF46yx2uk99b0JV6ARL/o+9x7D/caJwDlV8+8YUqoDnfYfHI6cxzGlsc6Icr7pvKqY0B9qnNmdqb8BzHuQahsZOr7Tg2eDyOLLk+QhoBi4EAOBOoTni2n/Y9hjeXDWTeWzdrJjoei8uhFhcZMCFAPDg7SuzIq49t+lIffuB8Wy0thudH6QQQjQgAy4EgEN3B/H4yLmmusYl8XMZX/AX1u7r1NCEEMKJDLgQKY+OJjnA67Ps8ZmY6Nn323+3bj/eWCUluhBi8ZABFyLl3KrmmOjnGnzB3ztHtDYhhOg0MuBCpIwW+jLlGV/ww+PX8cz1k0nBEa1NCCE6TVe5kQEWOUfCg1V5t4/LWKXAW2qTuxhLv8Bco5K2yDHMtYB0UiZ9GI/nCXflclNgbRHXEIBfW+Zexly/WDKTqidJAEvKwogbXD3Wlacy25oTLoxHqYE3BrAWb/duxfuqb8zZB3MXi4hrFEvqEZoEx3cMGxP77nI5flMxtzBLvj/Wh8+tzgR+rywJCN3f464V6v7FXJ2YK1fcghtZaD17trRy3u1yF2OuYr5jqsTtjbp4eVy5WFtVcgx77rCkJb5kWtz92OFG5vGg1AxciJSfHT7dFHHt5Ia9OL16JnGJb4YuhBCdRgZciJRDdwfx6K1UyJZGXPvBwOH69vWV5iXz9RX5ggshFg8ZcCE8NM6692d8wftrpcUYkhBCAJABF6KJ9eXsLHumPDta26N4c/U2CCHEYiADLkQD+yaGkw/pLHt/gy+4O1rb4dmNCCFEB+gyFbpbGcnUq75kH+NMzUgD2QcmFfBsa5vSnXRCxL8APKpyAhNS+jTBMbi61EV+jkQjWYz36rphiQuynO1v8AW3FifX78V7Bk8DAA7cHsZjN8/hxKYDM2OMTZNS3TlepkInXyBL6hGa/AQIV3xb4inQSh+5nPt7tTSZie/8wu4pRqh6O9kWlgiEqdDbqQRvRdHerr7Z/1Ko2tyXzIQnGgms91wP5inULrU52x/gdsbVls+jSTNwITxkleiPjl1KPkzP0Kdn7EII0WFkwIVo4KdvvDUr4tqpdTMhU8+tGkg+TEdrmy4LIUSHkQEXooH779zAB6+8khTSiGurKzMBXmb7gjeXhRCiU8iAC5FhMt+TfEhn2ecbZtlZlboQQiwWMuBCZBjLzKoby09morW9mHlHLoQQnaKrVOjWGkzVZg+5lVi0TEXN1OblQOW4D6YeZxpctj+DxVoHgJgESmfnUSbK4F6PErwnClMy9+bcsex7SYz7yKNWDlUZu9SuazOz7LWVqfp+B24P49GRizixaX8SrQ05PDdwGAduD3O1OblWlgQ5bkVtHgqLGc/iXPuIAq95O1X2obSmQg9TY4fGEGeqdSBcuR4a4956FM6sb6pCJ+fHFeL8vKmnDr3mpB1PLPRQb6B21fu2ueyPYqELEcDeOyPJh/Q/Z+9dv9I8O2MXQohOIAMuRIYLazYlH9LVilc27G7avrYy5S0LIUQnkAEXYhbNa1YvbtyDt9ZsqZezMdEbo7UJIUSnkAEXIsMTg2e8vuDnM77gL2dm6EII0QlkwIXIcHAOX/Cs8Ofkhr1NM3QhhOgEXaVCr1mDO5XCrHoWV3aKxCr2bZsk9WUirGZ958LF6cHq9FZgw8oT9WqeBFYvRnxUTAnL6iskJnIhdiuljadvHqTd3UfOuA8Yz/WmnSUx0c+t2oxq7W0YY/Hu62fw3NbDiBElM3QYvLZmJ+677V5Kt4Hqaqpmb6NK29bIfUvU6T6VNot5HpPvLyJ9LCa+84vJ/cnijoeq1mOPUjpUhd6u/QFPTHfyFGmvQjw07jiJw+5TupO2qoGqctYHsw0AV8D77kMXmoEL4cDnC37wzg383KXXksL0DL2s3OBCiM4iAy6Eg3UZX/BsebJQTD6k78HrynUhhOgQMuBCONh792byIV3/66uWm7b7ZuhCCNEJZMCFcHC30AvA1mfYf7L7KN5eu7m+fV1lomn/bFkIIRYaGXAhHBwZu4bIosmV7K+2Hqxvz0Zr66tUFmGUQoiVTFep0KvW4FbZEQudxfH2qtDd9SWiDqwQ4WyJtJP3/DRiCvVQjTFTrRPhuLdvdlCRxUInymOA/ypk9dVAxWnOG7+ZbXBX14g6fd/oCB4evoiXNu+daSI2dWXunfoMPQKsxX/cfRSPjVzEwTs3ZrXFVOg09jdTwbr3bgmmaGexsY3n7lyKMc9DaSUWOlNwM1U5VVa3MRY6VcDTMYXHQmf3SLAC3vN/zOOUE8V3C/kw2DHseRTah+/a8vjps+t8/y2agQtBODpyNflQj4k+Ut92ePS6d4YuhBALjQy4EISLmZjoFxuU5gdvD+GRmxczRyy9maUQYvkiAy4EYazY6y0fHbmcfJieod8ZgRBCdAoZcCEIa0uT3u0XVzfP0OULLoToJDLgQhDePXgWuemkJgBe3bSryZVsrChfcCHE4iEDLgThvttDOHozXSY3BjUT4a+33lffno3OJoQQnaS73Mhig8Gp2UOuEu0Qc/3ybWP1NRaRn1DqQL4G5i6WM9x9oUh+svXn3OfHkpb43G6KObfDU4EksWAJRZhbis/thrkoudwzAH4eORs7xzBa6Ee5lrievfPqOfzltkOoRcl4Xt60G2+u3jYrqUlExG3tdKdibYUmR2DJTFg9wF17mKuazyVtKRLqNlWtuV0TebIPjytXmxKH0P8lX0IRsi20D3Yf+BKptMtdzJfMJNRdjLkrM3cxXzKTCjmm6qj3mR7NwIXwsDYzy15bnkkreuD2sHeGLoQQC4kMuBAedt++lXxIfwbvvn3Tu39WqS6EEAuFDLgQHi6t3Zh8SF9LvLZpZ9P2xhk5AKzLlIUQYqGQARfCS/MLqFc278LZtQP1cn1Gns7Q99z1z9CFEKJdzNuAG2NyxpgXjTHfTcvHjDE/NMa8ZIw5YYx5p+OYw+n26b/bxpjPpNv+qTHmSsO2p9t3WkK0h3ddOwdjbT1kqjUGpzdsq2+/tLbZF/zVjTtdzQghRNsJUaH/KoBTANam5S8C+IK19o9T4/tFAO9pPMBa+yaAY0DyAwDAFQDfadjln1trvzzfAVQscGNqtiSPKYx9ynF2TEyOYRrcSuze0pdzK1FbITD/CFWaA0CBHJQjx1DFKe8iWPlMlbZMIduCwj9UCTutat01cgvvPXcK3z/wQHoxDIqlCqZqBQBANTPGlzftxhurtuLA2HC9LkeU/GxMiwlVlLN/GHAFPFOus/1buR6hSv7QxCStHBNTpfTCJxRh6nSuZg8/7+w9P7O/ux2PU4yn77Ax0UQxXqW7/38/C1OOs/3LnmvL2nLV+5K+zGsGbozZBeBDAJ5tqLaYMebrAFydo5n3AThjrb0wnz6FWCpMFYrJh/RJdHn6vTiAn7jqn6ELIcRCMd8l9K8A+Bya5z6fAfAlY8wlAF8G8Pk52vgYgG9k6n7FGPOKMeZ3jTEb5jkWITrKnYyyvLG8f+wm3nv+VFJIZ+irJWQTQnSAOQ24MebDAG5Ya1/IbPo0gM9aa3cD+CyA3/G0UQTwDIBvNVT/FoD7kCyxXwPwm+TYT6bv2E+UYkW+Ep1nTcYgZ8vZGfqlhhm6EEIsFPOZgT8B4BljzHkAvwfgKWPM1wH8IoBvp/t8C8AsEVsDHwRw0lo7OF1hrR201tastTGAf82Ot9Z+zVp73Fp7vCdSrGnReXbdTrOMpS/meivlpu3ZGfpt+YILITrAnAbcWvt5a+0ua+0+JMvg37fWfgLJO+8n092eAvCWp5mPI7N8bozZ3lD8KIDXAsYtRMcYL/YCsPUZ9nP73oFz62Yyj801QxdCiIXgXmKh/xKArxpj8gCmAHwSAIwxOwA8a619Oi33A/gAgE9ljv+iMeYYEjHcecf2WZTiGs6U7syqty3EVs6R3y4FcknyRAmYM+52JmrueOCAT1XO+nDXszH5YsD7FOouWFM+dSdV5waqPi2JKe2L9d6K4nUu9gwPIzoIxKkRjwH87bYD2DkyCgDYPjoK7EJ9hl4o1TBZLcyMl6jQWQz4TqjTWTxyGr/cc13ZMUwxzJTjrcSGD/UuYGMK9ZwAuKqcx/APj0dOVdfB/2NhY/X1Qes9ccddtHKfhz5DfGOqtkltTvdvQYVeDfSwCTLg1trnADyXfn4ewOOOfa4CeLqhPAFgVqJka+1/HTZUIRaHvaMjODJ4Fa9vm/Hxbnz0jBd7khoTAdbiB/uP4KGhK9g3qqAuQoiFQ5HYhJgHh4euJR/SWfau6RjpAO4bGUJkUXcli43BiR37Oj9IIcSKQgZciHlwdV3q5ZiuJV9ZO+P1uHf0Jo4MzhUGQQgh2osMuBDzwOcLDgDvGG6eoe9smKELIcRCIAMuxDxYUy55t9dn5I4ZuhBCLAT3okLvONZYVEx1Vn1E1IFRC79PanCrxy1ri8Vh9/TBlOss/i87C0Okwf15rn7MkU01on7sLbgHtbYw+3uYiwpTihK1OVMlF3xfa+A1ZGQvx9FLF/F3u/YhjpKWTm3ZjtNrNmPPaOIjPpqZkY8W+upK9DyJCZ4nKnQeQ3y+o58hIheEXVvm0WE8eQX4xXVvYAr4FpxJgmktFrq7vl1xylvx6AiNq87zDXhU6IHjDVXG+whVqIeq033bmBCcxmdnY6I98zwWoWgGLsQ82DM6gkM30mVyYxBHEV7atae+fXWpeYaeLQshRLuRAReiRe729NQ/bx9LfMJnorVVFmNIQogVhAy4EPPEN8uerPuCJ2tjz993CBfXKya6EGLhkAEXYp5kZ9n1MoB9N4dgGnzBrTF4sWGJXQgh2o0MuBDz5Nq69cmHdJb91paZvN970mhtQgjRKbpKhQ4AsUvbR1TdsFwHaGhEct6zC/amk8VaB7iil8X4ZmrzXiIp96m0e3Puvlfl3fW9RBHtU3cypW+OeQsQxWmuAzHBQ8gqXd/cug3n1g1g161EiX5gcAintu2sz9C3jN5GKY64apdIUXNEYexT5oaqzRksbjtrHwAs+d+ganOC7/xaUTK3C6b4Xug45f5jwmK605wGnljhPPZ3e9TprRB6P/ti3DPFPr+GYefhmx0zL6VQTxPNwIWYJ49cvpT8+JpeJofB+U0D9e3X161LPqT/hfWyEEIsADLgQsyTXbdG8FNvp1lzUyPeX57JDd6oSneVhRCinciACxFAqZCmCXXMslfJ91sI0UFkwIUIwDfLfuTyJURxXH8H/tbWbbi8Qa5kQoiFQQZciACys+zG8q5bI7h/8HpSSKO1vbJrdyeHJ4RYQciACxHA9rGx5MM8I67pPbgQYqHoMjcyS1xW3L9DYp/LgQ1MWkLr3Q4aBc9vozxxeytE7vr+HKknSUv63LlBAPAg+ixQP0tAEnsi9XtyqbQHj+tLqJtJKOPF4nRHgLX4m4P34+D169h1K0kf2p+ZofeXStTtpkqubZ6cA0uKAvjdvFwwd5UK6YIlXgEAQ86PuYWxvkPPoRVCXYGA8GQm7Uo0kvTdHtcsdq91IpkJq6/6Hs+B59fK/z2/tu7946Xl1QpAM3AhgtgzPNzsSmYMXt09s0y+bbR5hl4vCyFEm5EBFyKAXbdu4f7r1+j26+ubfcHPbN3aiWEJIVYgMuBCBHLf4I3kwzxm2W9t3YrLGzZ0YlhCiBWGDLgQgWRn2fUygKOXLs1aYr8wMOBqRggh7gkZcCECGc8oyxvLu27dwrveao7W1lcqQwgh2k1XqdBjWJTN7IdhniQtyXlOL4Jbqs2SkOTJ/kxtXjRcCt4TqDbvI7LuIvn51UMSUgAAyX+CIklywhS4Na+aN0x9TFtqIQmCaVPiBKagjq1B/1TzPWitaVLsThWalepX163H0fjSvPuwnr5Dx8vUuZ5bxEktNMsCfAlWyP4tJGvpBEwlHppoJDQBCcDV46Hq9NAEK95j2P5EZU+Tonj6ZorvUPW9755iLfHvKWx/r8o+oG/fna8ZuBCBPHTpclPEtTPbtuLKxpn33L4ZuhBCtAsZcCEC2TlyC/ddH0wKacS113bvqm/3RWsTQoh2IQMuRBtonGVvzURr65kjWpsQQrSCDLgQLeCbZU82RmsD8KP7DzYtsQshRDuQAReiBbKz7HoZwJ7hm7NcyRqX2IUQoh10lQqdETHluOVK8BxRlbdLbd6X432vojHM3fW9pKleohwveH6WFYn8mMb/JVpN3y8/tq3G5KvkgBCl5jRUdc26JvvPpdq9vm79dIeAtTizZSseOX8RALBj5BYOXruOt3Zsbzo221OoYjjnU9QGnndo7OhO/NL3KoYDzy8UT2j/4LjcTNlN1due9lnfoQp4pvj2K+Dd9aH3bYWq02nXFN43O8LzrCD1LF8Eg/XNzhvgKntntec6aQYuRCtk/qmySvT9N5qjtW0ZU0x0IUR7kQEXogUevHQ5s0wOXBrYVN8+uK45Wlu9LIQQbUIGXIgW2DFyC8ffOpMUUiPeW54J8CJfcCHEQiMDLkSLlAuF5INjlr1avt9CiAVGBlyIFvHNsh/MRGs7m3lHLoQQ90pXqdANDPJ29pAN+R3C1OkAkGeqcqJC74mIOp0EdmZKcwDoDVSbs9jmTG3u+1UWqtZkymdfbGqq6GWKV3JAjfbAicj3wWOCtx47vW+qlCmXUUljom8ZHsO+azdwdue2erS2V3fvxtaboy33x+JiA+Eq+8WklSuei9w3SegMJDSWto9QJX+oohxoIR55oBLcp0JnynXeVlg7rajQQxXtTO0NhKvNQ1Xrvr7ZeF31ioUuxAKwZTQs4pregwsh2okMuBAtMtVTTH4ep7P+Fw7dh6sb19e3z4rWNqX34kKI9iEDLkSL7Bq6CYPmiGuv75mJuJadocsXXAjRTmTAhWiRHSOj2H91kG6/sb7ZF/zc1i2dGJYQYoUgAy7EPbB/MBNxbZTPss9u24qrUqILIdpEl6nQ3cpypjbPeX6f9KLgrO/PuS9JD5EaMkW5T4VeJMNiqvIcqacqR4+6kqmuLdE6cuG4pxOmsCS7M0UmEZT7IW1FZLz3GmP7emPENWtxfd06vCO+DAA4dP4KXtu/JxmSMbDG4sKmAWwZXrpL6T7lLCNYzdvCNc+FxiMP2ttPqMKfjZRd2lZU6EyNzc6bqc2ZQjzp211fJX2z2On8vGnXwTHPWX3J48oSOi72fG7lMcX+Z6RCF6KDTPT20PL2kVEce+NsUkjfk/eVyhBCiHYgAy7EPdA/R8S1cjFd0UmXE25sWLvQQxJCrBBkwIW4B46cv9IUce3C9i241uBK1tL6mhBCzAMZcCHuge0jo9h7LRWypRHX3ti7s75988jt5MO0yO3W7U4PUQixTJEBF6LNNL4HH5peMk+X0M9v27wYQxJCLENkwIW4R7LvwZvKmSX0c9u3Ni+xCyFEi3SXG5k1KNjZ7l89xCVsteGxp3sj92+XfuKz1UfcwvrJFfR4kQW73dSIbwhzs/J6ApFj2JAMTcDg6yQUdx+tJOIITTjA/ORYshaX+9z6kTvAAdSXydeP3MFUmtTkwLlr+PH+PengEley8wMD2DA0Rr+/dv6q5m5F7euDeFJS2Hn7/mdCk5C08/5k9w5L9BOa5MR3bjx5R1hCEeb26UsoEuyqRurZ88uXnCfUXazC+vD4E1bIxWL3Tpm0lWvB3zXkf9z3v6oZuBD3yM3MMvnl7TPL5FtvjuKRjCtZr3KFCyHagAy4EG3mwvYtGNw0s0xeLqYrRKmBH96wbjGGJYRYZsiAC3GP3H/+ajK7TmfYMAbXtmysb5/oLTbtny0LIUQryIALcY9svTmKo2+cSwrWAgboaYi41j/VHH0tWxZCiFaQAReiDVQyEdduNkRcG7jVnFa0WK50dGxCiOVJV6nQcyaHDWbVrPqCcf8O6SFKcwDoJdLSVQV3PUtAwvSHPnVnqEC2xhJxkIbyPhUsOQ+W7IMlOfEli2CJE0JpJZkJ+8bD2wo7YLynZ1Z5WsU7XuxNvvQoSXjyypED2HF5GNtGbgX17PNeYKp5pvQNvQd9VyM0cUjEFMbe82tP363AxssGlQtsn6m9AV/ikLDvlSVMaUWFzpTdoUp3X98sCQnrmyVSqXkk3FzRHvbfwZPatKBOdxziu8c1AxeiDfRllsUby1sHR2Aw847cGoMz+3d0eohCiGWGDLgQbWBjJmRqvQxgy81R7LpyYzGGJYRYxsiAC9EGRjY2+4Jf2THQtH3n1eHkg8PACyFEK8iAC7EAXNqxBTcafMGzBr5eFkKIFpEBF6IN3Hduti/44NYZX/DJjO93tiyEEKF0lQo9bww2FmbHPWexmHs8st0eIhUtkJ80NM5vCzGXQw9hokgeg9qjfiSNESE/bcunjGTjZedNR8vilPtiZpN6piQOjVvNFPYbh8bwjtfP4dSDB+q+4PnJcl1F3DPZLHKLrV997MIXc5xdE6Y+bmcs9FCFPxur13miTXnVW4mRTq8t+ScgjizIUTk7hym7WUvsnmL3LVNvA1zxXQlUtLNrzuKX+/t2N8bjrfM+mEKd1bO+Gb4Y6SEz59jTr2bgQrSJSiZkauMy+YFzV2HiuG45r+7cjKEBZSUTQrSODLgQbcK3TL55eAw7rgwlBWNgowhn5UomhLgHZMCFaBM+X3AXeg8uhLgXZMCFaBNZX/BCJmTqLAM/qZjoQojWkQEXok2UeouJuih9B/7GA/sxNDCTOnRWsJdb8gUXQrTOvFXoxpgcgBMArlhrP2yMOQbgtwH0AqgC+GVr7Y8yxxwG8M2GqgMAft1a+5WGfX4NwJcAbLbWDvvGEAHod0hxcywmcQvKVZ8y0kUrYt5WlOsh+BS7TNVaJVJRX3xqRqjCmYk1WTstxUinyudwpTtj4PoIzMM2GbcxsADO7N+BDUNJMpObjb7g1uLy9gHse+tKQN++uM5hw5r55AAAIABJREFUAw6NLe6Nb9+m77sTs4l2xk6nHivk/NjDtpXcAdy7IExt7nsWsXjrTFXOlOOsD7Y/AJRIJ+z5XCUPC5+Cm6nNWVsV6+6c5YvwweKkR45/jnbFQv9VAKcayl8E8AVr7TEAv56Wm7DWvmmtPZbu8ziACQDfmd5ujNkN4AMALgaMQ4glycDwGLZfag6Z2vgsyD4XruzcguGGGboQQoQwLwNujNkF4EMAnm2otgCm/WTWAbg6RzPvA3DGWnuhoe6fA/gcWpvICrHk2HEtEzK1YZl8vyPYy42GYC9CCBHCfGfgX0FiaBtn858B8CVjzCUAXwbw+Tna+BiAb0wXjDHPIFmOf9l3kDHmk8aYE8aYE1N2cp7DFWJxuJUJmXqrwRd8YHgMR14/lxTSYC89JQnZhBCtMacBN8Z8GMANa+0LmU2fBvBZa+1uAJ8F8DueNooAngHwrbTcD+CfIFl692Kt/Zq19ri19niv6ZtrdyEWlamMa1i2nA32cksx0YUQLTIfEdsTAJ4xxjyNRLC21hjzdQAfQfJeHEgM87PkeAD4IICT1trBtHwfgP0AXjbJg2wXgJPGmHdaa6+Hn4YQS4NsyNQscxl4IYSYL3MacGvt55Eujxtj3gPg16y1nzDGnALwJIDnADwF4C1PMx9Hw/K5tfZVAFumy8aY8wCOz6VCtyCKPCbTa0HWGqoHDY7v7TmGwdpiKlgfNF4x66WN6gSqrg6MoewbE+uDnXfUxhOMYbDrzDWcO7gTNkq+nKs7N+PaxvXYNJwo0QuTzb7hxcnyLDUu86pgymMfeRbHmwQeZypmnxI21Kuincr/tsXdbwF6rUjM87iFe41dW6rsDtw/1HsB4NeW1TO1OVOaA+Fq8woJTM/2B4AaGXGNqM0r5L+gSjIw+L7viKnQHd9HTMaTtNM6vwTgN40xLwP4PwB8EgCMMTuMMd+b3ildLv8AgG/fQ19CdAWbhsew9XL6OzQNmXrxvpmQqetvNvuC18tCCBFIUDYya+1zSGbcsNY+j8Q1LLvPVQBPN5QnAGyao919IeMQoptoXCYf3dTsC3595wAOvD3bF1wIIeZCkdiEaDO9mZCp2XIj13ZtwU35ggshWkAGXIg241sm33Nmti/40Db5ggshwpEBF6LNNC2TAxjcOVDftml4DPe/dj4ppL7gxTmylgkhhAsZcCEWmGu7Njctk1d6mn3B6wZfCCECCBKxLQVcAeipC1QLmQtCE2t0gjx1K3LXt9NVphVCXblCXYG8BB/kHmygx1vTMTvfvorz9+9MKw1ggBtbN2LtjcSVbKqn2fd7srenKXEE8YhBjrh+AfyaV9g9EujK5XMVo99fYDKaVpL8sD7YM6EVt0HWliXuYohZopHgrukjjLl/hbqLtfK9Mrcwtn8rCWRYghCWnIS5izFXMaB97mI1tr9x7w8AMTnGOupjszBuZEIIBxuGxnAgs0xeaFgm78ksmc8V/EUIIVzIgAuxANQyIVNvNyyTr82I3Arl5uAuQggxH2TAhVgASpkQqY3lSm8xWYtPjfuZh/ZhZLNcyYQQYciAC7EAzFombyhvuj4CgxlXMmsMLjVEaxNCiPkgAy7EApBdJs81LJNvGBrD1otDTfsvpkhSCNGddJUK3QKoOQR5hiVmaKEPps5lqsh2kiPyXNZzqKq7UzDVPIOpWtt5Gu26JvNNuDHZky6TR0nI1HMP7cPmC0NYP5Qo0TdeHsb1vVvqlnvtzdv1RBRUbO7p3Ab+AlhMJXg7v9hQD4YaGVMr9wdLLsOSmUTkovtuKfaVhyr/Q5OiePsg+wd7mbTQt8sLCQAqRFHuU6GzZCNMbV5B1b0/UZvXyP7JMe5trmOYYh3QDFyIBWHjtdnL5FcOziyT384Ee7nREOxFCCHmgwy4EAvA+qExbMosk/sY3L0ZtyRkE0IEIAMuxAKxeTqtaMMy+TQ7354dE/2mYqILIQKQARdigbiTWSZv9AVf7wj2opjoQogQZMCFWCBKfUVvuZoJ9jKmmOhCiAC6S4VuLSoOFWKB6DhZHGjA98tlESXcTL1KYiu38+dXqArXq8YODH5cDVSvtgJXSocRcsnzE80zaovmc53IBHuZ6C2iVAN6cu72vIphFq+bHbAE3dZaUcCz27CdcfdpH1TJ796QY//GLSQvYDHdqfKf4PsfC72GXOnO4pf7+mbnFwZTmgNACe4IiDWiQmdq84pxr5x5VehkW83MHpPPA0ozcCEWiG1vXU38HtOH0fDuAYxtmRGq9WQMfLYshBA+ZMCFWCDW3RjDpkupkM0YIIpwrcGVbM1Ic7CXelkIIeaBDLgQHaTxPXhW5HZzl3zBhRDzRwZciAWkOOlZJs+82hravRmjW+QLLoSYHzLgQiwga4b5Mvn2M1cTRVCDL/gt+YILIeZJV6nQAbcCkmn9fOq9ULVmJ2Khg6rNiSKT7e/rop1qcwLXXrqpEGkpU6L6MG36XltRwLti2TeFTLUWQzsHsOXUFQDJO/K9r53HhYf3133BC6Uy7dt3Zmy4nY4Bfy+wmPhAuGg+VCntu9VY3HimKmfqdLa/D3bdWd6E0PNu5f4IVZuH7g+04B1CrkdEPAIAICL/UaVAtXkF7noW7xwAakQB70IqdCEWi8wz4uaezbjdsEye9QW/I19wIcQ8kQEXYgHZcnr2Mvno9pll8rmCvQghBEMGXIgFZO2NMex65XxSmF4mbwiZKl9wIUSryIALscDUepqXye8OzCyTZ33B85X5vxsTQqxsZMCFWGDKmWXxxnKlp5gotFLjfuHoPqUVFULMi65SoVu41ajVFtTK7QoGHbfQN1NMsjFViaw1b9z70+ZbgKubeSd5sonFPq4SKWyoEhXgMaJDaaVvS+6FwmRpVrlmLaLYYO3VEeBx22TEL9+3A2tvjM0ek+fU2KZ23QqtqJXbpRxvpY8a+QLZPehT+hpyr7O1EqYcL5DpUivfUaiinarvW+g79JHHlfG8IRoLPbDeB1OhF23BPSbmX8O8ETxK86otOetdx1jP00gzcCEWmP6h5mXyqDTzT7ruxhg2XRhq2n8J5hkRQixBZMCFWGBqfc3L5Fcf3t/kSlaPl54a+LU3FRNdCDE3MuBCLDBrr44AaHAliwxu3D+T1OTOQHNMdPmCCyHmgwy4EAvM2htj2HjhRlNdo04hK3KTL7gQYj7IgAvRAdZfbF4mr78Xx+yEJ0IIMR9kwIXoAOOZZfLxBl/wbW9dTWTTqXEf3j2AUbmSCSHmoOvcyCrx/B18TBv9qZiLEBuNz60hzxIRsID8zE2hTQk3fDCXsBxJsAIAVeJ2w1xDKoFJEFqBnXVoko5WXJ0KMM6QqdP3zpobY9h0aRg3920BjIGNIly+bwdWDc52JWsXzA2pnQp4dq3Yv4bvPzt0psHaYu5izO1s+qh2wFwvW0kUw74/1lbodwGEu9yxa0gTBvGu6fOzRtpqJdlURO8qdzITtn+OmNEc3O5oSQ9uF7OadbmRKZmJEItKIXCZPPteXAghssiAC9EBBk43L5OPZrKSZQ283osLIeZCBlyIDrBmcAzrL6YBW4wBchEGG1zJ1gw3B3tZLV9wIcQcyIALsUg0LpPfzYjcRnYNLMaQhBBdhAy4EB0iu0zeWM7KVIb3bMbYFinRhRCc7lOhBwStNy0lOeF9u4hbUD+GqnBzge37fpWx0TJVK1OUkzwq3l6YqpzVM4+DVhKNMPW9R0wf3DfzLphW3/dlYqL3D42hksqDB05fwfXDO5MvzxjAADe3bUT/9Rkluk+t7EsM4cIQXX6oihngngrsWoWOFWC6YA47v1YIVWPTvltwYWBOIyxNRqiXSSvqbXYa3CPHXd/K7LFm3b1U2TPHc+eEPrtN4IgjG36Gxsw+xveNagYuRIeYyCyT39o9s0y+ZnAMO145lxSsBQyQl5BNCOFBBlyIReLW3i24s3VmmbxWTP1GHcFehBAiiwy4EB1i0+mryXridFITY3B7+8b69qzvt3zBhRA+ZMCF6BCrB0ex9eXMMvnUzDJ51ve7OKUldCEERwZciA7iWyZflfEFz5WYVEkIIbpNhW4tSna2qpD9CmmnEpXBYvP6YOONyRYeQ9ld34r2PlS96otTzlSnTH1cbWOceeotQI5h+/PY6eH3VON9WHHERJ++luXeYjKgyADW4srD+7D2/A2smUdMdBrrmpwh8zpohVpgW210DvH1ErS3z7uAxvIOjRXO1Ni+6xd4rWiM9LBmvISqzUPb8cGf6ewZwgdVIyOIWxpZGDnjjpMeWZfPEb9BNAMXooNkleWNvuBrro4AaHhHHhkMH9oBIYRwIQMuRAfpz/qCD8+ETF0zOIb1528sxrCEEF2IDLgQHWRyc7Mv+Nju5pCp6y4NJx8cBl4IIRqRAReig2TfyI1mfMGzwV4m5AsuhCDIgAvRQTa9OdsX/M6OGV/wrMgtWxZCiGm6S4UOiwqqs+pzNFp4++SuEVECMpWjL8ZwzNoicX5ZNN8c+f1VcMTTnYbpGZlC1gQqan0wFTqrZ8pxn/Lfp1B3ERrPOt/Cb97GGOn9g6PY/NI5DD12oO4LbibK9ZjouYlmkZu1qG/zEaq+j9uoQmfObl519QLTQthxSmhTzKuC4YvHz7we2vW9tvM74qp8N76Y+BXyLKwQhbjLLgBcaQ4AVeN+stZIW2XjjstQgbu+ZsLdQPPomVXn86bSDFyIDhP3NPuC19+LI52h1+L6O/CxvZsxvnV9x8cohFj6yIAL0WF8y+SrBkex7sJQUjAGyEUYOSxXMiHEbGTAhegw2bzgc2Ud03twIYQLGXAhOkw2L3g2ZGrWoCutqBDChQy4EB2m1peGTE3fgd84tr/pPfesYC9D8gUXQsymq1TogFv1HRPVoA+mKo/IbxqmfGZxc30xeFnfXG1OVPZEqcnaT7a5Cb2Cvl9+VDkb2EeoshrwfU+hWuIwxTwAFMhViTIn3nt5GHjnwRkjHgFDh7ajeP0WAODuwJpkR5PERB/dM4B1r1/yjtY3rnbs78Ow+PAs9nfbeua0M5q1JdeKnTfbP1SdDvD/JZ5PgbRD6r05DQLbYrB7rRSzJx5QJk/DMnlSlYhCnCnKAaBCjqmSYyqm5O7DutXmrB4AIvJMjx3n7fNo0gxciA7Tf30Ua84Nznv/O/u2YEJKdCFEBhlwIRaB1dNK83R20tewTL7BEexlfOdGVzNCiBWMDLgQi8DU5jR8qsMXvH9wFAMvnksKabCXnIRsQogMMuBCLALV/qK3XOtJ5SkOAy+EEIAMuBCLQn7CP6Oey8ALIcS8DbgxJmeMedEY8920fMwY80NjzEvGmBPGmHc6jjmcbp/+u22M+Uy67X8zxryS1v+pMUbhpsSKYf2bV5pCpt7Z2yxUyxr4uQy+EGLlEeJG9qsATgGYXsv7IoAvWGv/2BjzdFp+T+MB1to3ARwDkh8AAK4A+E66+UvW2v813fY/Avh1AP/YNwALtzuQT2bPYA4MJtABpZW+WTITliiD9cEcJHKeMeVIohN21swFpI25InhihpYSxYQdE+oGmKeJc3hSiorj+y5cu4VVF25g/MC2esjUm4e3o3BtBABQvDGWDjxtdKqCSszvTepyR7YwV8OWcluQ74+5QHHnIU7oUiG7UqHJbnwwdzF2zVnf7Dolfbjr6fMrNJGKp29KoPtcpYWkRCWSIqdEXLlYopGSmaJ91EgfzP2rat19s/rY8js9MsSNzHGMJe7CwDz/L4wxuwB8CMCzje1ixpivA3B1jmbeB+CMtfZCMijbGJ1iFdprE4ToOqr9M5mIssFebj12AJPb5EomhJhhvjPwrwD4HIA1DXWfAfAnxpgvI/kh8O452vgYgG80Vhhj/ncA/xDAGID3znMsQiwL8uMlWu6/chOAbQr2MnZkJ/quj3Z2kEKIJcucM3BjzIcB3LDWvpDZ9GkAn7XW7gbwWQC/42mjCOAZAN9qrLfW/pP0+H8L4FfIsZ9M37GfqIIvhwjRbfRmQqb2NviC910fxaqzmWAvWqMSQjQwnyX0JwA8Y4w5D+D3ADxljPk6gF8E8O10n28BmCVia+CDAE5aa1n4qX8H4L9wbbDWfs1ae9xaezyP3nkMV4juYGpL+gYqXSYf37u5aXs22EuvYqILIRqY04Bbaz9vrd1lrd2HZBn8+9baTyB55/1kuttTAN7yNPNxzF4+v7+h+AyANwLGLUT3k5lR392/tek9d9bA18tCCIF7S2bySwC+aozJA5gC8EkASN3BnrXWPp2W+wF8AMCnMsf/hjHmMBLB6AXMoUAHAGssDVrvgiUm6QSRDVd38iQkbhUiS3LiS2YSqjanql2ferWNSl8XPj8Brth3K0LDVes+BTxR+JNDVr1xGWMP7kl+RhsDGIs7OzcinyrRK309TftX+npQIYrUUG8I5glBE9F4vlN2DL11SDu+MwhVrocnrwmnRr6L0EQqrSSWYU82E/i9+vr2PUdcsHukSq5TxZNopGpIMhPy/K+AJTPhCUWYejxYnR67X+2y9kOxnjsqyIBba58D8Fz6+XkAjzv2uQrg6YbyBIBNjv2cS+ZCrBR6r9/CupNnMHb8YD3ueWPI1NxE8wMgWxZCrGwUiU2IRSTuKSQf0hlSaTpGOoCeoWZf8KjEZxNCiJWHDLgQi0itv4eWs77gY4/dh6ltGzo5PCHEEkYGXIhFxLdM3nd52hc8TSsaGdw5sqvDIxRCLFVkwIVYRHzL5L3Xb6E/6wsuhBAp96JC7zgWFjWHcpGpzWvBelBO1XDFpIu859LmyLaKdffB1Oas3qcdZfG0Q2Np+4Sz7brqoTHSk21hsc1rTJ1uWoivTzwPqmRMReRR7ssnFz8ygLUYe+w+9Jy7hp7rtwAAPReuY+K+bfULnh+6hYqjPaZCp6pkeg7h581iYIPE3Wfq5laU4+weaSVHQXDfpJ6pq0Pj0ifb2LPNTY7lUyDXid0f6UY3gZeWxTxn/xfJMe5ryOorLEa6nfD04daUsNjjTIXOVOIxeZ4DQBzP357ccyx0IcTC0HN5GNll8vEju+vby5tTv/D0PXi9LIRY8ciAC7GI9Fy/hZ6z1+h2n8hNCLGykQEXYpHpu3Aj+ZAucxaHZhKWyPdbCMGQARdikalM+347lslXvXEJqMV14z61bytKciUTQkAGXIhFx7dM3nP9FnrPp0p0Y4Bc1PSOXAixcuk6FXrFEQvXdOB3iC8erROfFJyoOGOi2mX7V0gnFcuvR570weIV85jg4bHeGSxOM8OnMA6ObU7U5kzN7iMmp0HV1dPVE1OZatukjs+OvdJf8MaRng954sHQGuy+ZTGwg1rxwr6lVhTfjFBvFu7xQNTNnrZywTHu3X0wNbtPhc7yOeTIMa3kDwiFnQev5/c5i1XO6mMaC53s71Ga1+JJZ711KNet5w7RDFyIRabvjQtNy+SlfdtQ3raxvj2XMfB6Ly6EAGTAhVh0CtdHUDyfKtHTZfKJI3tmtmeCveQbRG5CiJWLDLgQS5DG9+CVjC94ae+2xRiSEGKJIQMuxBIgmhUTvXHZvPkdYnl/8xK7EGJlIgMuxBKgMJSETp1eJq8vmwPof+MSEDdEazMG5V2bF2OYQoglRJep0GOUzGz1nmmrojaMyKP4phDhJ2vLGreytErUiUw9CgA91v2VV4OVom1UlpLxtqZKJnGXjftasdjKoWpeH1TF3HDapS0NvuDWYnLvFvT8+GzS5/Vh9J98ExPHj8wY8cnJpnbbOd5Q2tV3K54NoTHPW4nTz86PxfJmY2L3oJfguONumJrd543AYpiz74mdd5nGNefXg3kW8XqWF4Lfg4Z45JjAZzqPhc78LfgxPsW5C83AhVgKZJ595f3bUWlYJo97ismH9D14dbOCuQix0pEBF2IJ0PPmee8yeZwJ9hL393Z4hEKIpYYMuBBLgPz1m+g5+UZSSI14NDkjbMuK3LJlIcTKQwZciCWC7S0kH9Jl8krDMnlW5GZKMuBCrHRkwIVYImSXxRuXzeO+npnldQCTjx1pekcuhFh5yIALsUSIxptDpjYukxcvDyUfpo14ZDB1ZG8nhyeEWGJ0nRtZhQSabxeRCXRHMAVndUwSOSSQIPeGfB3ElSRHxsTcWwAgYi4PtD7cXSw0YQRznGBuSL5EI8yVhbmLVU1YUhBf3+Euhc19m6GbAA7Wl8ltqVR3tYmuD6Fw9ioq9+2s728bxsOSskQscU7g9+pz12Lfn6HuZcxtMNyNLPQ8WrmfQ93C2D0SnBAJoC6nwUlLiKsmc/FK2gr7nthzh7mLsXvWR564wQJuUafPlZE966dobp7/v71zj7Hjuuv49zdz777sOqka22vHya7tpE7SpE0d98WjtCkFGtLwElIpUCQEBf6iRQIUFSSeEn0gtaUSKCpCVLRFFFSoqNLSUoJCIaRJnLRJHCfxY2N7vWsbaGLv69478+OPmfue39k9N3fvnbn7/UhXuzNnZs6Zc8/M786Z7+/3M/pcst3FLDe15GBWQVYddj/xCZyQvNAxTb52+CZUp1/VKC7PpfHSUwMf1t+LE0K2JDTghOSE8Fya97tlmnztptlGeVQXtaUGPqIvOCFbGhpwQnJCaeESwpNnzPJukRt9wQnZytCAE5IjSo5p8qAjLzghZGtDA05Ijoh3pa5hGdPk48+cBqK4Ydyrs3va3pETQrYWBVOhK2rYuAo9RLZC3EVsCP4CT4GsS/0YGt1e1rHM9SXNVktadfRTzWsm4nAey0+9aiuD/VTM7jZZClk/FboTy4vAIKtN0VT7GKhtG2+oncPFSyjNzaN2YF9i4MMAqzfPIFy85PAisL4Lf4V/v1Czbvt7tZTgvu21lM+9nLelKrc8G3pJ+lIzxqeZ+KiP9wTrWJbXgW/CIJcHSL+uS1ffljCeuX7SUI+HgZ89UYcnklWW5b2kDpvHJ3BCcoR4TpPzPTghWxcacEJyROn4ibZp8mhmL2rT1zTKu4K9LPG9OCFbFRpwQnJEuHAR4dx8spBOk1du2t8oL3XERC/RF5yQLQsNOCE5p3WavLar3Re8OrNnGE0ihOQAGnBCckawtGIvd2iwqvv3okYlOiFbksKp0CPtjjsbWvHITb2kC2MfSyhtiDstpTkAjGu28MhSm5eRvT401pccitOaqeb1i/ccGLGVXfTybfTrOL6q614Uw1aMezOOvsWliwBubEyT49JFVKWCkpZQOn4CldccTH56iwAiqFx7DSYWLmYeylZd+8Xx7ieBMdZcY6pf6nFL3RxL/87brMP0hLAxx62hlLZihVeNeN0uXPewzDpQMerOXl8x1gOA9ml8ujyRLCV/CZOZ68uGar0cZt/PK8F2s+5aKVtZvhZd6Vq3Upkzj8MncEJyRrwzfaKu+4LPNBOYlBYuofzY08lCPeTqCnODE7IVoQEnJOfEs9chmt7ZXDGR+oqnBr4R/IUQsqWgASckZ5SOn0giCtWfsEUQXzvdKGdMdEIIQANOSO4IFy4iPPpkstCYJm/6eztFboSQLQMNOCF5ZLxjmnxnU2kedPiCY81foEQIKT6FUqELADHUl1lYsYqTY/n9drGOFSH75ulSS1pq8+2GylEMVXlkKcdNyTxQNRSya4Yi1DpvSwXrwq3g9jhOD2pl81hGm3zjWSftsuJQ+513hO6Y6NHUGNYkeQqvTYWJ8Q4CQBXVw7dgdW4OoaFE92lrP7HO24qv38uYshTf/n3ev5j4VvzyqBcluGarqM3zMy4NM2+CYxxInzwSLLV5RZfNfcz7jqcPikuFXpbs+22ghtoc2bkqzBwWxvEBoCbZwtOydCvg16rnzOPwCZyQHNIZE711OTjXknJUBAgEtUMHB9k8QkgOoAEnJIfIxfRpOp0mbywDCBYuQE7ZvqGEkK0BDTghOUR3pglM0nfgOnNdW3kwdzbdMDHwwcX/GVjbCCH5gAackFzS/jJTZ2cQT+9qLncEe2kVuRFCtgY04ITkkOD4s0Acd/iCNxOX6FS72EXpC07IlqNQKnQLS5kYGLHCAVvlaKnTrTpsdbqtaq0YatSaoWa0fmWtGQp4S2kO2LHNfVW4rl9+dh/2El884zgOpXm/Ynn3Eos58lQAWwrjWGJgcR44+jhwx+GGEY9WriBOFb26vNTVrjVsPDd4SYZ36VvjoxcvBd9x6zs21yvL3N6Iq27GSFeHstpyuDDGjhWPXz1jqgP+ngpWzHPLU8flJVQ1xnKshtdBD+PZsg9WzHNrfJoqdCMuPQBUJXufSoZy3WXH+AROSF4ZTy/mdJocO1vCqR4/DkRR0xd85nro9O7Bto8QMlRowAnJKx3T5K3LsrgIzKVKdBEgDKGHXj3AxhFChg0NOCF5ZXnFvdxB53txQshoQwNOSF7p8AXHWkf0Jk8DTwgZLWjACckrk5NNFToAvP526O6W99xdwV4uDbiBhJBhQgNOSF45l8ZAboRMDYBDh5rldVGbEeyFEDLajIQbWS/YrizZkv1Qst02wh660HK3WLPcETyTQlQdSROsJB3eOHOG+LndWC4u1nfkcj+xzs9yibG+PzHGgTrcliL1S1YhyE5o0BiDiy9ATp2AHLgBgqTLFRE03U8QJQl+0jLMzqAyfTWweH7duiPP83ZhudxZeCfocOA7nq1x4Eo04utGZrqWeo4PV93+31I2rmvJSkJiuUdZbbXO23IvA4BanH1tuNqbiZ33yByH1pgaj7NjLVjj1rqvucgen/ZJ8AmckByjc6eTv+k0uV680Cw7fgyIY6gqJA32ItfuG0YzCSFDgAackBwjO5PwqZJOk9eXAQCL5xEffRRAauBFoCsUshGyVaABJyTPTE05lyUN9pJp4AkhIw0NOCF5ZnnZXb6OgSeEjC404ITkGD1+DIiixjtwmdkP7G4mNeky8OsZfELIyFAwFbpkqsR9FeUAUMoIGg8AIQyFrKFeLRld6FKnW8HvIzNhSrbk20yK0oPS3DvXW3akAAAU30lEQVQJiEPdaR3LFdw/C0tx6jo/Uwnr2SdWW11qZQuX2jaT1q974QTKc88jPHAIIgINQ+DQIUQLLwAAgovnEeK1DQMfrS0h0tWWQ/kl7THHv9jXkpWMwzpvqw7L08OFrypZjOQdruNYKmrvBEfGcaw2AYAY953A2Cc27nm+iUkA+95mXX/mdWl4rFj9l+yS3Yext6eHfd5qJJ0xExlZSWqMBCuu+1TN8GapSrf63jU2+QROSNFonSZPg73U34GXXv8myPS1Q2oYIWSQ0IATknc6UodqyzS5njuT/K0b8SBAeOjWgTaPEDIcaMAJyTnxxQUATV9wXFxslOnCPKJTzw2jWYSQIUMDTkjOCXZOA2hxFZs50L7B3CkATQNfN/iEkNFmwwZcREIROSoi/5wu3y4iD4nI4yLyiIi8MWOfQ2l5/fOSiHwgLfuoiDwjIt8WkS+KyNX9Oy1CRpdg9iBkem9jWXYmCU7qBr5u8Akho42PLPjXARwDsCNd/giA31fV+0XkrnT5ba07qOpxALcDyQ8AAOcAfDEt/hqAe1W1JiIfBnAvgN92NUAQYEy6/VxtFbr9+6Sslgo9u0ssVbK1vavuQB0S7gxq4qd27QUzvrGhvOwF61i+CnjX9qZy3RHDPAvbu8BWSluqVl+Vdtyhao2feQxTt9wODYJEiS6CeO80quefBwCMT423j7apycaxVf361lLZB+offdtSXVv91EuscMvTxFYf+8VtB2y19Ea/v8b2lgrdoRAPAqPfDWW35S1gKd3d9ylrHz96uU+51ONZWOr0yHEca6zHyLYNVWTHhrc8ciLHPce6F/r21YZ6SUT2AfhRAJ9uq6tpzK8CML/OYd4B4ISqzgGAqv6LamOkPwSAQZwJySBeOIPKY98E0BoytUXItrRk7UoIGWE2+jPn4wB+C+1ppj4A4KMicgbAx5A8Qbt4D4DPG2W/COD+DbaFkC2HTCSZkBrT5LuawVxqx59oC/YSztyIYJq/hwkZddY14CJyN4ALqvpoR9GvAfigql4H4IMA/tJxjDEA9wD4QkbZhwDUAHzW2Pf96Tv2R3yd+AkZFWRqe8fytsb/8cJZRHOJEl1EgDBE6dBrB9o+Qsjg2cg78O8FcE/6nnsCwA4R+RsA70byXhxIDPOnjf0B4F0AHlPVxdaVIvILAO4G8A5t+Mi0o6r3AbgPAMbCHc4s1ISMKrp0pWPZPW3eauAJIaPJuk/gqnqvqu5T1Vkk0+DfUNWfQ/LO+wfSze4E4HJG/Rl0TJ+LyI8gEa3do6oM4EyIg+jieQAtecHXVtvKdfmKc5kQMnq8nFjovwzgEyJSArAK4P0AICJ7AXxaVe9Kl6cAvBPAr3Ts/ykA4wC+lr7Xe0hVf9VVYYAQ2/SqDTfQrULPjkc+pv7xmLMQV7BwXywVcw9V2MpSK7ayoRh2KSzNfYw29TEcgVW3rfo04h73oJy1PBIiQyltxnvO6NtgcioJmRoEUFWMHf4erJ18ohETXS6cQQlHGga+cuEMqrF/bnA79ra/Ct1S39txrm2FuBUn3VJXu7wFfOmX2rwaZz+nWHkZAKCm3bGxAZjKZ/P7M8agM/68VYdxLO+Y/30kNrwtAocXjfU9WedRM/qqlzjz5j6e93QvA66qDwB4IP3/PwDckbHNPIC7WpaXAbwqY7sb/JpKyNaldu40xqCNkKkaBBi7+Q6spAa8tDPxCxcRqCrKM69G9alvDbPJhJBNhpHYCCkA8cIZVE8+3bbO9WO9vP8mhNPXb26jCCFDhQackIJQPX0cQPM9eO3CuUZZ5ZmjQBw3k5pIgNK+/UNpJyFkMNCAE1IQwl1JmtC6L3h9GQCihRew+tiDAJrBXuIVakMJGWVowAkpCN2+4O3LwXh7sJf6e3FCyGjyclToA6ekIa6Ou1XoUU+KYUutmf1mMTSVxNmu6bUeYi5bdZvrDSVj2aGkjyW7vZaq3FR1O/q8KoZy1iDsk/Lfha/aXD1jpwNuFbUPVhxoXb7cuaah6I7iGjDVkSdgagq12O+7sAgMtbcL6zysgEyuONDWPoGhDLbirVuqdReW2rwWr2avN5Tj5nEcdQfGtWEppa3zs7Z3xZ+3c0xk3wutmPG214F95vY+nt4kjvNTY+yYfWWp041+KiPb0wmwcy2sZt5D7PAnfAInpCBUjj3aFjJ1bPYWlKZnGuVxh4HvXCaEjBY04IQUhNrCHCqnEyV6PWTq2M1NT87ahbMAWvKCr/n7gRNCigMNOCEFpvU9eDC5PQn2kr4Dnzr8dpT2zFi7EkIKDg04IQXCFTK1evZ5oCXYC4IA4zcfGXALCSGDggackAIRpb7f9WnyqMUXvHZ+DpWTTw2lXYSQwUMDTkiB6PQFL8/e1FZeOX0MQGuwl7MDbB0hZJAUyo0slABXyUTX+qrhWhA55PcWvTjLZK/1zzRiuYtZxMb5Wa5tAFCxXDeMqgPDyaXqcvkxXbP8Eo1YSU6sZAr9xHIJcyVssM7PTJrg6U4lCLrGVXn/zShNz2J1PonSFnbERA937kFkJcTwIHJcSpa70SAIxEgoYrgyWm0VR7IWNcaC5S5Wi7LdyyxXJ1f3Wd+dNUasJC79cnEE/BPb+Ca1Sfaxrg2/87DczgD7+gsk2yxa3/ckslP3jhkJswCgZozbckZiG9f1xSdwQgrE2rFHu0KmlvcdbJQHU69o275zmRAyOtCAE1IgkpCp/w6gNWTqUqOcvt+EbB1owAkpGDI+mfzNiIm+euzhtmAv4/tvQ3kPk5oQMorQgBNSMFwx0WvnT2Pt1JPJ+jTYy8TNbxxk8wghA4IGnJCi0eEL3rXcAd+DEzKaFEqFHgCYCrsVkJU4W0Jd0x5U6JJ9rNg4lllDD6pP/5Qs2bgU8P5Kdz/lOOBI5mAoSC3lrBpqV3H0lKVc91XG+64HHKpkz6QX6lDOAsDa4imM3fbmpqvY6uW2xBDR0ott20dL30XUkdTEV8nvq/510YsXgZWkw1b+Zx/HUhj34LCCKK5kro89Vf/qSOZjJaMJjOuvBs+6+/i99hPrGnB5aBgFfSM0kp+IkVTKlfDJSh7le23wCZyQgtEVMvWOd7S9565dOAOg6QteufDC4BtJCNl0aMAJKRiVs8+hM2TqxC1vapSXdl8HoCXYy67rh9FMQsgmQwNOSMGonT+NtRPfaV/ZMgUcTO1oK+pcJoSMBjTghBSQelrRZsjUM42yeOmlobSJEDJYaMAJKSClXe3T5PVpcwBYPfbfbb7gkwdeh/KeA4NvJCFkUymUCl2RLRYNDeW4GOuB5pPLRvex6oiM48SO30aRobC09uiXOt2FSzGZhUu9asULr8TLmestdaelTnfFBrb2sbDaaqldLeU44Ip5nq04jY31nYrxOq0KVZns8AWf3N6IH702/yxWTj6ByRsPJzHRwxCTt7wFlfnnG9uroYLt7Te979jpoW5DJW6p0y01rzq8J3yx4nVb36uFGMfpBZeXROb2Do+H9bwhuus2Yp57HseF9X2HQXcMcQAIPO8Hzrpz+LybvxYRQtYlXvabJud7cEJGDxpwQgrIytMPOafJtcPAxx2+4YSQ4kMDTkgBqZ4/idWTiRK9HjJ18pa3NMrrvt91A1+lLzghIwcNOCEjQus0ed33u67pmNh/21DaRAjZPGjACSkone/BW6fNOyWX4wdfh7E9B0EIGR0Kp0Kvxt2Kxl40jmZsc2O9pULvJ9Z5xL0EajawYqGXNFu1WzVUn76qdcBWyK5G2YKswKg7MFTrgK1o7xeW0hxwqZINdbqhNo+NGNudCtzK4glM4a2NMbu2eKqx79JTD2Lqtu+HBmHyFC6C8nU3YHX+mNl+wFbaihVD3EHf4lm79jHiUFtYKuZeMM/P8FQQK+a/Q7VuXWV23/rFNu9NIW7kR7DOo4e2inEdB4bJMuPJO27b1ngLDA8Gy8NlUqcy17vukTXDEyPIHM+u3BaEkEJS3j0LoHWa/LWNssr887j8rfsBpD9KJUC8cnngbSSEbB404IQUlY7ZoomDt2Ns7w2N5WAieTJoxERPDT4hZDSgASekoKw8/Z9AHDWSmogIxq+7qVHeHRP9qkE3kRCyidCAE1JQKudPdE2TRytXGuWdMdHpC07IaEEDTkiBCSa2AWhOk4/tnmmUVS7MAWgKM3UtO5QtIaSYFEuFroqVuFu5qIZKu+ZQb1tqbItQ/ba32gTY7YoMVWbNUDO66rCwlJGxWG3KVkv6xlwGbIV4LNkq7Wq8kl23Lpl12LGxLXV1/37DWircODb60NjeWh9r9/iQyW0dy9sRpyroYGIS0BgShFBVvOINd2H5+YdRmX8Otr7ZX8VvjgVv9XG2+reXus06rJjqrrqNMaKWR4Kl7DZuIW4luJ+S39y+j/HIzXFrfq9+4yM5lrXe+i42P2NEtkLcvqeO65h9LEsBnzHWrG2TMkJIYYmXXzSXV888BQCNd+QIQmx7zVsH2j5CyOZBA05IgaksngLQnCaPW6bJK/PPYfm5R4bSLkLI5kMDTkiBCSd3JNPk6TvwHW+4G2N7b2yUr556HEDTwNcNPiGk+NCAE1Jg1psmH9u9H0CryG3/4BtJCNkUaMAJKTDrTZN3+n7TF5yQ0YEGnJCC45om7xS5EUJGh0K5kUVQXEF3ogDT/UrsYPmWNN9K6hEa25d76EKrvVVj/ZpkJ0eoieFu5HDxstxPrH2sOlxYLlslGc+u2/ieYjHapNnuZQBQiyzXs+zzsN3Ievlt268kD9nuSZrhRgYA5dT3W0Sgqijvnmm4OF158hvYdtvbgTAZp5MHD6O8Zz8q8896nYM7SYa/m1Dm5sb37cb6nnxd21zj3KjDPJblBOXn2pYU+Sd+8TmOcx9H4p7sHSz3suwTFKcrr5/7nDU+XePWt08iw93VcmFzuSpPGi5mQYa7cta6Zh2EkEITTl1tLlfmn8XKiUcBJAZewhK23fq2QTaPELJJ0IATUnCi5e96bd9p8AkhxYQGnJCCs/TkA9Co1ngHPnnwDoztfXWjvNPA+xp8Qkg+oQEnpOCsN01eWTgJoEXkli4TQooNDTghI0jrNPnY9AEATV/wyQOHh9ImQkh/KZYKXSK8GHRP//WioC5p9qmXJVsdWDZUg5GhQCw7kiPEhirTTDRirLcSjdQcilpLSWnVYan1XSrY0EiIERnq6lKQrU7vRTlrEcVG38aVvtVhqcf91crWedty5Wj5/zKW68dp32/yhiMY23sjKvPHN16DuhLn+LbXuDYMlb2b7H0shbO3QjzZyRMra4mh0nap740EGr6JX7wTjfSA73i2EpMAZg/adVvX2ACSnPSCaQMykkq58mjxCZyQEcA1Tb705ANAHDWitYkIJq6/dRjNJIT0ERpwQkaAzmny+jIAVOaP48WH/xFAauAlQLT80uAbSQjpKzTghIwALl9wAAgmkrzhWQaeEFJMaMAJGQHWcw1bz8ATQooHDTghI0C3L/gRjO091CinLzgho4eoU12aL0TkIoC5YbejB64BcGnYjSgI7KuN0dVPpVfuPRiMNx+t45XLF2svLr4AAMHU1deUdlwzUy+Lll9ciF66eG5grR0uHFMbg/20MQbdTzOqujOroFAGvKiIyCOqemTY7SgC7KuNwX7aOOyrjcF+2hh56idOoRNCCCEFhAacEEIIKSA04IPhvmE3oECwrzYG+2njsK82BvtpY+Smn/gOnBBCCCkgfAInhBBCCggNuCci8tMi8pSIxCJypGX9O0XkURH5Tvr3zpayPxaRMyJyxXHcsoj8dbr/MRG5t6XsARE5LiKPp59dm3eG/WEY/dSyzZdE5Mn+n9XmMKQx9RUReSKt9y9ExC87xhAYdD+JyJSIfFlEnknr/ZPNPcP+MKTxtO7+eWRIfXVHuv55EfmkiPjmbmmiqvx4fADcDOAQgAcAHGlZ/3oAe9P/bwVwrqXszQD2ALjiOO57Afxt+v8UgNMAZtPltrqK8BlGP6XrfhLA5wA8Oew+yHNfAdiR/hUA/wDgPcPuh7z1U/r/29P1YwAeBPCuYfdD3vppo/vn8TOkvnoYwFvSa+/+lzOmCpVONA+o6jGgGVO6Zf3RlsWnAEyIyLiqrqnqQ1n7dB4awDYRKQGYBFABUNiME8PoJxHZDuA3ALwfwN/150w2n2H0larWx1YJiXHKvRhm0P2kqssA/i2toyIijwHY16fT2TSGNJ42sn/uGHRficgeJD+e/ys9xmcA/DgSQ+4Np9A3h58CcFRV1zz2+XsASwDOA3gBwMdU9X9byv9Kkunz331ZUy75ot/99IcA/hTAcl9bmQ/6PqZE5KsALgC4nG47CmzGtQcRuRrAuwH8a78aOmQ2pZ9GlH721bUAzrZsdzZd1xN8As9ARL4OYDqj6EOq+k/r7PsaAB8G8EOe1b4RQARgL4BXAnhQRL6uqicB/KyqnhORVyCZ7vx5AJ/xPH7fyVM/AdgB4AZV/aCIzHoec9PJU1+lYwqq+sMiMgHgswDuBPA1z+P3nTz2U/oU9XkAn6yvGzZ57Ke8kqe+QjJt3knPs1804Bmo6g/2sp+I7APwRQDvU9UTnru/F8BXVLUK4IKIfBPAEQAnVfVc2q7LIvI5JINj6AY8Z/30KgB3iMhpJON6l4g8oKpv66WN/SZnfdW44arqqoh8CcCPIQcGPKf9dB+A51T14720bTPIaT/lkpz11YNofw2zD8B8L+0DOIXeN9Ipti8DuFdVv9nDIV4AcKckbEMilHhGREoick1aRxnA3QAKo7DuZLP6SVX/XFX3quosgO8D8GxejHevbOKY2p6+i6s/Xd4F4Jl+tXvQbFY/pcf+IwBXAfhAv9o7LDazn0aNTbxPnQdwWUTenL4KfR8A5yyAk17Vb1v1A+AnkLy3WAOwCOCr6frfQfLO4/GWz6607CPpPnH69/fS9fcA+IP0/+0AvoBEMPE0gN9M128D8CiAb6dlnwAQDrsf8tZPHXXPolgq9EGPqd0AvtUypv4MQGnY/ZDDftqHZHrzWMtxf2nY/ZC3fnLtn/fPkPrqCJKHsBMAPoU0oFovH0ZiI4QQQgoIp9AJIYSQAkIDTgghhBQQGnBCCCGkgNCAE0IIIQWEBpwQQggpIDTghBBCSAGhASeEEEIKCA04IYQQUkD+H6qigosJHsd4AAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 576x720 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "f,ax = plt.subplots(figsize=(8,10))\n",
    "ATL06_gdf.plot(ax=ax,marker ='.')\n",
    "show(dem, ax=ax, cmap='inferno');"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Plot data"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 45,
   "metadata": {},
   "outputs": [],
   "source": [
    "ATL06_gdf_sub['srtm_height'] = refdem_sample"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 46,
   "metadata": {},
   "outputs": [],
   "source": [
    "ATL06_gdf_sub['diff'] = ATL06_gdf_sub['srtm_height'] - ATL06_gdf_sub['h_li']"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 48,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "''"
      ]
     },
     "execution_count": 48,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAX0AAAEICAYAAACzliQjAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjAsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+17YcXAAAgAElEQVR4nO3dfZxU1Zng8d/TLzRvDTTvDY28KCBvCUKDqAQdR0dwZyPqjsLMRB0yQVyzGyaZVTRjktmJbpvddWfcSQIYiePEDpKFoMlAEs0oKtCNXYAidAiv1XR4x1JbQRq6n/2jbmF11a3u6u7qqnurnu/n0x+qz71V95SeevrUc849R1QVY4wxuSEv0xUwxhiTPhb0jTEmh1jQN8aYHGJB3xhjcogFfWOMySEW9I0xJodY0DfGeJqIqIhc0cHnfiwiY7r6On5iQd9HROSwiNzUjvNvEJH6qN93i8gNXVI5k9NE5Dsi8pNM1yOWqvZW1YOdfR0RuU9E3kpFnTKtINMVMOmjqpMyXQeTfUTE4oiPWE/fGJOQiDwsIn8QkQYR2Ssif+z06v+fiPxERD4ClgCPAnc76ZR3nOe+LiLfFZEtTvkvRGSAiLwgIh+JyNsiMirJqtwkIvtEJCQi3xcRiarjIhGpdY79WkRGRh27lLJxrv2LqGt/16X3HncdEZkALAeucd7HBx3+D+oBFvT9Z6qIvCsiH4rIiyLSPdkntjc9ZHKbiIwHvgrMUNVi4BbgsHP4NuD/Af2AZ4EngBeddMrno15mAfAlYDhwObAV+DHQH6gFvp1kdf4UmAF8HrjLqQsiMp/wH5w7gEHAm8BPE7zG94FPgKHAvc5Pm9dR1VrCf9i2Ou+vX5J19iQL+v5zFzAXGA18Drgvo7Ux2awJKAImikihqh5W1QPOsa2qul5Vm1X1XCuv8WNVPaCqHwIbgQOq+qqqXgR+BlyVZF0qVPUDVa0DXgOmOuX3A/9DVWud13yCcMdoZPSTRSQfuBP4tqqeVdU9wL+04zpZw4K+/zytqkdV9X3gF2RhozTeoKr7gaXAd4CTIrJaRIY5h48k+TInoh6fc/m9d5Kvczzq8dmo540E/klEPnDSLu8DQvibRbRBhMcwo+vt9h4SXSdrWND3n6xvlMY7VLVSVWcTDq4KPBk5FHtqWiv2mSPA/araL+qnh6puiTnvFHARKIsqG9GO62TNcsQW9I0xrkRkvIjcKCJFwKeEe+ZNCU4/AYwSkXTHlOXAIyIyCUBE+orIn8WepKpNwDrgOyLSU0SuBO5px3VOAGUi0i0Vlc4kC/rGmESKgArgNOFvmIMJD5q6+Znz7xkR2Z6GugGgqj8n/O1jtTOT6D1gXoLTvwr0Jfxe/pXwgO/5JC/178Bu4LiInO5UpTNMbBMV/xCRw8Bfq+qrzu/fAa5Q1b9McP4NwE9Utczt+cbkMhF5Ehiqqm6zeLKWBX1jTE5wUjrdgF2Ep2VuINwJWp/RiqVZm+kdERkhIq85Nz/sFpGvOeVTRaRKRHaKSI2IzIx6ziMist+5meOWqPLpIrLLOfZ09A0WxqSbtW1vEJEvODc9xf2k+FLFhPP6nwBrgP8NvJTia3ifqrb6A5QC05zHxcDvgYnAb4B5TvmtwOvO44nAO4TzgaOBA0C+c2wbcA3hKVUbI8+3n47/EM6xfuzyszHTdfP6j7Vt+8nFnzZ7+qp6TFW3O48bCN9FN5zwFKY+zml9gaPO49uA1ap6XlUPAfuBmSJSCvRR1a2qqsDzwPy2rm9ap6pPaPguwdifRINZxmFt2+Sidi2U5KyTcRVQTfimjV+LyP8inCa61jltOFAV9bR6p+yC8zi23O06i4HFAL169Zp+5ZVXtqeaxiQtEAicVtVB6Wjb1q5NugQCgY8I3zU9N/ZY0kFfRHoDa4GlqvqRiHwX+BtVXSsidxFef+Mmwl9vY2kr5fGFqiuBlQDl5eVaU1OTbDWNaRcRCaarbVu7NukiIvvcAj4kGfRFpJDwh+IFVV3nFN8LfM15/DPgR87jelre6VZG+OtxPS3vhouUG5NJgrVt41GBYIi12+sR4I5pZUwfWdLp10xm9o4Q7unUqupTUYeOAtc7j28E9jmPXwYWiEiRiIwGxgLbVPUY0CAis5zXvIdcHDk3KVdZXceXnq2msrquXc8Lp98ZibVt40GBYIiFK7dSWV3HC9V1LHymikAw1OnXTaanfx3hpVF3ichOp+xR4CuEFzoqIHyL9mIAVd0tImuAPYTXunhQw7dAAzwAPAf0IDzDYWOn34HJaZXVdTz6810AvLnvNHVnPmHZrROSeu7mzZsBBgA3Wts2XlN18AwXmj7LEl642EzVwTOd7u23GfRV9S3cc5YA0xM853HgcZfyGmByeypoTGtefLtl7375Gwe5edLQpD4Ys2fPBgioarnLYWvbJqNmjRlAYb7Q6AT+woI8Zo0Z0OnXtbV3jK8VFcQ34eWbDricaYw/BIIh/mz5Fu5bVU23/DyKuxcwc1QJP/3KrPTk9I3xqkAwxM4j8TvX/fvvTqYk92lMugWCIe5avoW3D4doON/Ex41NNHx6kW2HQ+w93pCSa1jQN75VdfAMF5vjZ/02NytVB89koEbGdE7VwTM0JVgO7VsvvUfFhtpOX8OCvvGtkp7dcIn5qHPMGL+ZNWYA+QlGUC82K8vfOMjS1Ts6dQ0L+sa3dh/9MOGx91o5ZoxXTR9Zwpol1zK8X/eE56zfebRT6UsL+saXAsEQL9Yk3qZ1/4nU5D+NSbfpI0t4euG0hD1+oFPpSwv6xpfWba/nYqLkJ/D24ZAN5hrfivT4JwwtjjuWJ3Rq6qYFfeNLseE+9uuwAhUbOz/oZUy6RaZsfuVf3qa4ewFP3D6lRfAvyO9c2Lagb3zpzmllRE/RP/VxI2Uxgd96+8Zvoqdsvn/2AtsOh3jspfeYNrLk0h2yjRebeXjtux1u2xb0jS9NH1nC3TMuu/T7xaZmJg7rG3fek9bbNz7iNmWzqVk52XC+xbfb/Sc/ZsHKrR0K/Bb0jW9NigryzQo3jB8c19vfZr194yNuUzbz84TBxUVx515o6tj9KBb0jW9FT9nMA0JnG/nPfzQ27ry12+vjyozxosgA7oxRJfTvWcjMUSWsuf8a7ohJZwIU5kuHBnTbtXOWMV4RCIb4WdSUzQJnMarpI0tYv6OebYc/693bDuXGT6aPLOFnS66NK//vt03hqVf38uHZC/TrUcjf3Dy+Q2vxWE/f+FL0EgwC/Kfpn20w8fC8CRQ635EL84U7ppUlehljPKViQy2znniVu5ZvaZGWDARDfOulXZxuaORCk3Lq40Ye/fmuuPOSYUHf+FL0EgwKTI4ZxHU2SKGpWVO2UJUxXaliQy3L3zjI8Y/Os+1wiLtWfDZQG+7kxD9n2+FQuwd0LegbXwqdbbyUthHn94i12+svfUCaFR576T0bzDWe96vdx1v83hS1cOCsMQPicvoR7R3QtaBvfKnh3IVLU9hiF1g73XC+xblNzWqDucbz5k4a2uL36Dtvp48s4cX73e/Qbe+ArgV94zuBYIgfvXXo0u+xPf1BLtPbXny7znr7xtemjyxh49I5zJ86jIK88B+FQb278fdfnNyuAV0L+sZ3qg6eoSlqTeX8vJY9nTumlZGf13LOTlOzTd003hab3mnW+IXVKqvrWL/zKBebw8dPfdzIt19uX/rSgr7xnX0nGlrcnfjXs0e36OlMH1nCP9w2OW6qpk3dNF4Wm96J7cwAbHzvWNzzUp7TF5ERIvKaiNSKyG4R+VrUsf8iInud8u9FlT8iIvudY7dElU8XkV3OsadFxD6Hpl0iPZ1oxT0K484bP7SYmM5+izt4AY4cOQIwztq28YJlt05gyZwxDO1TdOmmrNi0zbzJpXHPa29OP5mbsy4C31DV7SJSDARE5BVgCHAb8DlVPS8igwFEZCKwAJgEDANeFZFxqtoE/BBYDFQBG4C5wMaka2ty3qrNh1r8LgmWmXVbwyR2Y5WCggKAelWdaG3beMGyWyew7NYJCY//+dXh9aZWvXWQcxebmVTah/uvv7xdOf02g76qHgOOOY8bRKQWGA58BahQ1fPOsZPOU24DVjvlh0RkPzBTRA4DfVR1K4CIPA/Mxz4YJkmBYIiDJz9uUVY+ssS1wUfWMIkO/LFd79LSUoCzYG3b+MefX33ZpeDfEe3K6YvIKOAqoBoYB3xBRKpFZJOIzHBOGw5Eb2lU75QNdx7HlrtdZ7GI1IhIzalTp9pTRZPFqg6eaZHLzwOWzXPvFU0fWcI/zJ9Cfp4gQLeCvFbvzE1H27Z2bRKprK6j/PFXmPDYxk7vgduWpNfeEZHewFpgqap+JCIFQAkwC5gBrBGRMbiPl2kr5fGFqiuBlQDl5eWJt0cyOaWkZ7cWDWbxnDGtfq3986svY/zQYqoOnrm0Lo+bdLVta9fGTWV1HY/+fNel3yNjVv+44KouuV5SQV9ECgl/KF5Q1XVOcT2wTsP3u28TkWZgoFM+IurpZcBRp7zMpdyYpETuwlXCvXy3AdxY0xOkf6II1rZNBrnNyHn99133TTCZ2TsCPAvUqupTUYfWAzc654wDugGngZeBBSJSJCKjgbHANmdsoEFEZjmveQ/wUkrfjclqs8YMoKgwj3yBboV5ndonFC6tzzMSa9smg9xm5NwwblCXXS+Znv51wJeAXSKy0yl7FFgFrBKR94BG4F6nZ7RbRNYAewjP/HnQmd0A8ADwHNCD8CCXDXSZpE0fWcILfz2rzXRNsjZv3gwwALjR2rbJlMig7FOv7uWTTy9yy6ShXZbaAZDIaoReVV5erjU1NZmuhslSIhJQ1fJ0X9fatelKrbVr20TFGGMyqLK6jlVvHQQRFl03ulPTMZNhQd8YYzIkduZO5HFXBn5be8cYYzLEbeaOW1kqWdA3vlFZXceXnq2msrou01UxJiXcZu64laWSpXeML0S2kgN4c99poGu/AhvT1QLBELuPfnjp3pN0sZ6+8bxAMMTKNw+2KOvqr8DGdKVAMMTClVt5obouLuBbesfkvKqDZ2iO+WRMKu2TmcoYkwJrt9fTGLsMrMPSOybnRe9/G5HMEgzGeFFldR2/fCd+lY7eRfk8eutEm7JpzPodLbc5zE+whr4xXhUIhqjYWMu7Rz7gfIIefjoCPljQNx4UCIaoOniGhnMX+OnbdXx47mKL4388YUinl2AwJl0CwRB3Ld8St6lPtIHF3Rg/tDgt9bGgbzylYkMtK988GJfDjxDg/usvT2udjOkMt13cYp1uaGTByq2sXhy/RWKq2UCu8YzK6jqWv5E44APc38Ya+sZ4SSAYcs3fu2nvBucdZT194wmV1XX8/S92t3nezZOGpqE2xnTe0tU7Lm2Ikoz2bnDeUdbTNxlXsaGWR3++i/MXm1s9TyAtPSFjOqtiQ227Av7MUSVpSe2A9fRNBlVW1/H91/bxhw8+dT3ev1c3PvikkcifgsKCzm+cYkxXq6yu40dvHWrXc64fPzhtaUsL+iYjYlcXdBP6pJHCgjyuHzeIwcVF3DGtzPL5xrPa6sQkkq60ToQFfZMRydxqrsCFi81MHdGPB//oiq6vlDEdlEwnBkCchXYicxWuGNybJ+/8XFo7M5bTNxnRcO5CUueJ3YhlfCCZTkyf7gU8Pn/KpX2euxfmpT3gg/X0TQZUVtexs/7DpM5d/AWbomm8LRAMUXfmbJvnLZs3gT+/+jLGDy1O2T7PHWFB36RVIBji79a3/jW4qCCPEf17pmXrOGM6I5m7bfv1LOShW6681JanjyzJaEemzfSOiIwQkddEpFZEdovI12KO/62IqIgMjCp7RET2i8heEbklqny6iOxyjj0tIpLat2O8bsWmA63efJWfB5VfmcWrX7++ywP+kSNHAMZZ2zYd1drdtj0K83ji9ins/NafeKrzkkxP/yLwDVXdLiLFQEBEXlHVPSIyArgZuLSVkYhMBBYAk4BhwKsiMk5Vm4AfAouBKmADMBfYmNJ3ZDyrsrqO3+w5EVc+Y1QJjRebGdKnO/dff3n6pq4VFADUq+pEa9umI2aNGUC+4Br4h/Tp7qlgH9Fm0FfVY8Ax53GDiNQCw4E9wP8BHgJeinrKbcBqVT0PHBKR/cBMETkM9FHVrQAi8jwwH/tg5ITWZjfcflVZRj4cpaWlAGfB2rbpmFd2H6dvj0I+vdDE2Qstby6c69G7x9uV0xeRUcBVQLWIfBH4g6q+E/NNdjjh3k5EvVN2wXkcW+52ncWEe01cdpn3/lKa9mkt4AsQOtuY3gq51SMNbdvadXaJ3sITIE9g/JBiQmcbmT91OMtunZDB2iWWdNAXkd7AWmAp4ZTPN4E/cTvVpUxbKY8vVF0JrAQoLy9P5/aRJsXaGrgtSPONKW7S1batXWeXX+0+3uL3ZoWGTy/w/b+Y7ukZZ0nN0xeRQsIfihdUdR1wOTAaeMf5alsGbBeRoYR7OSOinl4GHHXKy1zKTRZbt72+1YHbkf17ZvoDIljbNh3Qz2X3tvoPPuWuFVsJBEMZqFFykpm9I8CzQK2qPgWgqrtUdbCqjlLVUYQb/TRVPQ68DCwQkSIRGQ2MBbY5YwMNIjLLec17aJkvNVmore7s6EG901IPN6oKMBJr26adAsEQe080uB5rak7PEskdlUxP/zrgS8CNIrLT+bk10cmquhtYQ3gw7FfAg87sBoAHgB8B+4ED2EBX1rtzWlnCY3kCSzK4IcrmzZsBBmBt27RDIBji7hVbOHfBfVXY/LzMpyxbk8zsnbdwz1lGnzMq5vfHgcddzqsBJrevisbvYqe09S7K59rLB6Z1eqab2bNnAwRUtTzROda2Tayqg2eIXQU8X6BnUQEThhbz8LwJmU5ZtsruyDVdqurgmRYpnpsnDuGZexLGWGM8b9aYARTk0SLw/8P8KZ6ck+/GFlwzXaqkZzfyRMgjvMBUJtM5xqTC3uMN9Oru3/6yf2tuPC8QDPGdX+zmYrMiAvddM8rTX3uNaU3Fhlqe33o47iYsCK+yaT19k/NWbDpAo/MdWBWeefOgp6eyGZNI5EYst4APMG9yaZpr1HEW9E2XCARDvFrbcp2dJrU9bo0/xd6IFW3+1GG+6eWDpXdMigSCIdZur2f/iQbe/6SRD89diLspq5sH7r41JlmV1XW8+HYdH5y9wNEP4tfLH1zcjaU3jfdVwAcL+iYFAsEQd6/cysVWFhW/YlAvnvxPn7ecvvGF1taL6lmYR2FBHheboe7MJ2muWedZ0DedtmLTgVYDvoAFfOMrrW1/ePZCMzi5/ciCa15dXM2N5fRNpwSCIV5xWSM/2uWDelnAN77SnoHZ1vL9XmRB33TKik0Hklpfx2btGL+orK5j1VsHKSpMLjx6dd38RCzomw5zm6GTJ/DE7VO4eeIQBhd3A+CVPSdY+EyVBX7jeZFc/v5Tn3A+wfTMiII8YcmcMb5K7YAFfdMJbssmF+XnsWrzIc5faOJkw2ebozRebGbd9nqM8bIX365r+yTHf79tsu8CPthArumEUw3n48rOXWxm/8mP2X/y47hjtmuI8bLK6jr2JVguOVq/noU8dMuVvpuqGWFB33RIZXVdXGqnNfl5rS+zbEwmtTZFM5of0zmxLOibdgsEQ3xz/S40ya57n+4F/PivZtoMHuNZrU3RBJg5qsTzSyYny3L6pt2e3FibdMAHWJYlHxaTnQLBEHVn4u+4jXb9+MFZ04atp2/apbK6jm2HE8/CeeL2KYTONtJw7gK7j33EvMmlvs19muwXCIa4a8UWmlqZqFOYZcuHWNA3SQsEQ/zd+sR5zyVzxliAN76ybnt9qwE/m9I6ERb0TdKe3FgbN0UzYv7UYb4f4DK5JRAM8esEd9MWd8/nub+6OquCfYQFfZOUig21rmmd4SU9ePCGK6yHb3wlEAxx1/ItJFoy6pF5E7My4EMSA7kiMkJEXhORWhHZLSJfc8r/p4j8TkTeFZGfi0i/qOc8IiL7RWSviNwSVT5dRHY5x54WkVY3XDfeUFldd2lhqVg3jBvk24B/5MgRgHHWtnNP1cEzrgG/f69uPHG7f/a77YhkZu9cBL6hqhOAWcCDIjIReAWYrKqfA34PPALgHFsATALmAj8QkXzntX4ILAbGOj9zU/heTBeITM9MxM83XBUUFADUW9vOPYluwvrbP/Hf+vjt1WbQV9VjqrrdedwA1ALDVfU3qnrROa0KiNx5cxuwWlXPq+ohYD8wU0RKgT6qulVVFXgemJ/i92NSbMWmAwmnZwr+vuGqtLQU4CxY284lFRtqWb/zqOux0NlG1/Js0q55+iIyCrgKqI45tAjY6DweDhyJOlbvlA13HseWu11nsYjUiEjNqVOn2lNFk2InPvo04bH754zJmrxnOtq2tWtvSLQUcn5edk3NTCTpoC8ivYG1wFJV/Siq/JuEU0AvRIpcnq6tlMcXqq5U1XJVLR80aFCyVTRd4O4Z7l91s+F29Ih0tW1r197Qr0dhXNnQPkWsuf+arOnEtCap2TsiUkj4Q/GCqq6LKr8X+FPgj52vtRDu5YyIenoZcNQpL3MpNx4VCIYInW1k/tRhvLHvNHkC0y4r4f7rL8+mD4dgbTunfHDuQlzZ2CHF2dSmW9Vm0HdmITwL1KrqU1Hlc4GHgetVNfoe5peBShF5ChhGeFBrm6o2iUiDiMwi/BX6HuD/pu6tmFQKBEMsfKaKxouf3bnSrSAvqwK+E8tHAm9a284NFRtqOfZhfMqyPTtl+V0yPf3rgC8Bu0Rkp1P2KPA0UAS84sxOq1LVJaq6W0TWAHsIfzV+UFWbnOc9ADwH9CCcJ92I8aR12+tbBHz4bE38bAn6mzdvBhgA3GhtO/tVbKiNm3rcr0cBD82dkPUzdqK1GfRV9S3cc5YbWnnO48DjLuU1wOT2VNBkxvYEu1z5eYpmrNmzZwMEVLU85pC17SxTsaGWZ96Mv9dkSlm/nAr4YHfkGhcVG2qpPR4/jzlP/D1F0+Qmtx5+RC6ldSJsaWXTQiAYYoXLByRP4Lvzp2RNasfkjkRTNHOVBX3TwopNB+JSODNHlfCzJdfm3Ndgkx3mThqa8Fhbm6dkIwv6pgW3m7GybWlZk1uW3TqB+VOHuR6z9I7JeW43Y1UdPJOBmhiTOmOHFJMXNR2lT/eCrF9YLREbyDWXRG7GmlrWl531H14qL+nZLYO1MqbjKqvrWPXWQc58cv7SXhCF+ZLTezZb0DdAOOD/xY/CN2NFb5Qi5MYiVCb7VFbX8ejP41eIvdCk7D3ekLNB39I7BgincGIDPoAIObEIlck+rQ3S5uIAboQFfQOEA3tBfnxzGDOwV872iIy/tTZIm4sDuBEW9A0A00eWMLWsb1z5otljMlAbY7rOkjljcnIAN8KCvgHC+c+3Y/bAnTGqJKc/HMbfEqVwil2WVs4lFvQNgWCIx156r8VNWXkCy+Zlx3r5Jje5pXAK83Njo5TW2OwdE94kOmYEt3xkieXyjW9VVtfx/df20atbPgX5Qu/uhUwq7ZNVS4N3lAV9Q4PLphJXDCnOQE2M6Ty3qZoLZ1yWNTu9dZaldwy/fLflJk9+3/Dc5Da3XL4tuvYZC/o5rrK6jvoPWq63Uz7KUjvGvwb0ir+D3G1f3FxlQT/Hrdp8KK7MBnCNn535JP4Ocrd9cXOVBf0cFgiG2H/y4xZlE4bmzgbRJvtUVtfx7h8+iCtvbXnlXGMDuTlsxaYDcWXTLOAbn3IbwC0qyOOvrh1lg7hRrKefw3bUxe+De4cN4BofCgRDPP5ve+LKiwryLODHaDPoi8gIEXlNRGpFZLeIfM0p7y8ir4jIPuffkqjnPCIi+0Vkr4jcElU+XUR2OceeFhG3DddNGgSCIU593DL3ObC4W06ldo4cOQIwztq2vwWCIf5s+RY+aWyKO6Yauw+cSaanfxH4hqpOAGYBD4rIRGAZ8FtVHQv81vkd59gCYBIwF/iBiOQ7r/VDYDEw1vmZm8L3YtrhyY21cWXTRuROwAcoKCgAqLe27W/rttfHrQ4b8RdXj0xvZXygzaCvqsdUdbvzuAGoBYYDtwH/4pz2L8B85/FtwGpVPa+qh4D9wEwRKQX6qOpWDf/5fT7qOSaNKqvr2HY4PrVz//WXZ6A2mVNaWgpwFqxt+1VldR2r365zPTZjVImldly0K6cvIqOAq4BqYIiqHoPwHwZgsHPacOBI1NPqnbLhzuPYcrfrLBaRGhGpOXXqVHuqaJLwvV/9Lq7sikG5vYRyOtq2tevUigzcNjXHH8vPE5t6nEDSQV9EegNrgaWq+lFrp7qUaSvl8YWqK1W1XFXLBw0alGwVTRKWrt7hOmc5l5dQTlfbtnadOoFgiO+8/J7rsT7dC1hz/zU53YlpTVJBX0QKCX8oXlDVdU7xCedrLc6/J53yemBE1NPLgKNOeZlLuUmTpat3sH5n/H/ykf175vISyoK1bV8JBEPc+cMtNDa5J/KXzZtgAb8VyczeEeBZoFZVn4o69DJwr/P4XuClqPIFIlIkIqMJD2ptc74mN4jILOc174l6juliN//v110DPsBTd09Nc228wZnZMRJr274RCIZ4eO27rsfyBJ64fUoud2CSkszNWdcBXwJ2ichOp+xRoAJYIyJfBuqAPwNQ1d0isgbYQ3jmz4OqGplL9QDwHNAD2Oj8mC42u+K3cevrRCyZMyZne0WbN28GGADcaG3b+wLBEHet2OKawwf47nwL+MkQr89jLS8v15qamkxXw5cCwRAP/iTA8YbzrsenlvVl/Vdnp7lW3iIiAVUtT/d1rV233zd/vosXquNn6uQL/IMF/BZaa9e2DEOWiuQ9E7GAb/wmUfd08vC+FvDbwYJ+lnps/a6Ex5bMGWPzl42vVFbX8W/vuo9J3T3DAn57WNDPQoFgiD3HGlyP2UCX8Ru3hdQAiovyeeTWidae28kWXMtCVQfPuJZbwDd+5LYTFkB+fp615w6woJ+F3Pa8vXniEPuAGF+aN7nUtfyGcXaDW0dYeifLBIIhnnmr5W5Y+QJLcmxdHZM9Ip2Vp17dy/sfN5KfJ/yHKaX844KrMlwzf7Kgn2XWba+nKWrJQXGms+XqXHzjb4FgiOWbDrDn6If069GNr9803r6xdpIF/ZrV42AAABOwSURBVCzz1r7TLX6/ckixfUiMLwWCIe5esYWLUTdjRQZ0rU13nOX0s0hldR3B98+2KOtWYP+LjT9VHTzTIuBHrNp8KL7QJM0iQhZxm+Vgc5iNX80aM8D9gMdXEfA6C/pZJHbWztQyu1PR+Nf0kSUsmRO/5HcuLwOeCpbTzxKV1XXsrP+wRVlxj8IM1caYzlu6ege/ePcoAnQvyGN4/54sum60dWQ6yYJ+lnjRZcu4RPObjfG62L0fzl1sZljf7hbwU8DSO1niw5jUztA+RfYBMb71+u/jt5Pcdvj9DNQk+1jQzwKV1XUcPtNy1k7vIvsSZ/xp6eodfHA2/q7ymaP6Z6A22ceCvs8FgiG+6bIY1ehBvTNQG2M6J9GWnnPGDuT5L1+dgRplHwv6Prdi0wHXdcZt2QXjN5XVdbz8TnzAL8gTC/gpZDkAn9t2KH5FzRmjSmzZBeMriZZPBpg8rE+aa5PdrKfvYxUbavng3MW48mXzbIMU4y/ff22fa7nt8JZ6FvR9KhAMsfyNg3Hl1ss3fnTuQvx6C1cM7m0Bvwu0GfRFZJWInBSR96LKpopIlYjsFJEaEZkZdewREdkvIntF5Jao8ukisss59rSISOrfTu74uwRfha2Xn7xFixYBfN7admYFgiEuH9QrrnzRdaMzUJvsl0xP/zlgbkzZ94C/V9WpwLec3xGRicACYJLznB+ISL7znB8Ci4Gxzk/sa5okVWyopfZ4/HaIS+aMsV5+O9x3330AsXkFa9tpFAiGuGvFFt4+HALCAWlwcTfb5a0LtRn0VfUNIPauCAUioyt9gciQ+23AalU9r6qHgP3ATBEpBfqo6lZVVeB5YH4q3kCuqdhQ65rW6dejwDY7b6c5c+YAxA6KWNtOoyc31tIUldlpBk5/3Mj4ocUZq1O26+jsnaXAr0XkfxH+w3GtUz4cqIo6r94pu+A8ji13JSKLCfecuOwy+2sfkSiPD/DQXAv4KdJlbdvadUuV1XUEgqG48mYNL6ts31q7RkcHch8A/kZVRwB/AzzrlLvlMrWVclequlJVy1W1fNAg2wcz4rH17nn8OWMH2lfh1Omytm3t+jORKZpNLv+l8qSVZZVNp3U06N8LrHMe/wyIDHbVAyOizisj/PW43nkcW26SVFldx55j8Xn8+VOH2Y0rqWVtOw3c9n4AGFjcjZ8tudZ6+V2oo0H/KHC98/hGPhsMexlYICJFIjKa8KDWNlU9BjSIyCxnZsM9wEudqHdOCQRDrjeujOzf0zaHTj1r22kwqdT9hquv3zTeAn4XazOnLyI/BW4ABopIPfBt4CvAP4lIAfApTp5SVXeLyBpgD+EBsgdVtcl5qQcIzwTqAWx0fkwSvvpCwLX8urED01yT7LJw4UKAKwGxtp0+gWCIV2tPuB4LnW1Mc21yj6jHtx4rLy/XmpqaTFcjYxLdnp4n2NfgFBCRgKqWp/u6udquI1M0m1z2vi3MF1YvvsbadAq01q7tjlwPq6yuc11Bs7go3wK+8aV12+tdA35Zv+4W8NPEFlzzqEAwxN+t3xU3DWRocRFV37wpI3UypjMCwRCv7z3pemzCsL4W8NPEgr5HVR08Q7NL5u2/3jQu/ZUxppMCwRB3r9jCRZdePsDg4qL0ViiHWXrHo0p6dnMtt4Eu4zeBYIiH176bMOB3K8jjjmll7gdNyllP36N2H/0wriw/T+ymFeMrbfXwJwwt5ru3T7HUThpZT9+DAsEQldV1ceXTL+tnHw7jK09urE0Y8AEGFhdZm04zC/oe9I01O13v479iiC1CZfyjYkMt2w7Hr60Tbd7k0jTVxkRYeifDKqvr+P5r+widvUAecL6pmQsuC5LkCdxpeU/jExUbaln5pvvigBBePnnpTeNtzagMsKCfQUtX72D9zuSWafnufMt7Gu+r2FDLC9VBGs43tXre+KF9LOBniAX9DKmsrks64N88cYh9QIznJdrrwY2ldTLHgn6GJFplMFZ+Hiy5/vIuro0xnfer3cfbPMfSOplnA7kZMqCX+zz8aBNLi1lzvy23YPxh7qShbZ5jaZ3Ms55+BgSCIX7xbsuevghEr31388QhPHNP2tcBM6ZDKqvreLX2BPl5QpPbreQOS+tkngX9DFix6UCLD4bQMuDn54mldIxvJFoJNlq/noU8dMuV1sv3AEvvpJnbWuKx/aIbrxxsKR3jG8mMT00Z3tcCvkdY0E+zddvrXRdSixBs4Nb4SzIpG0vreIeld9Js34n4fW6jXT6ol/XyjW8sXb2DX76beOqxpXW8x4J+Gi1dvaPN29IXzR6TptoY0znJ3FxoaR3vsfROmiTzAZk/dZh9QIxvvP77U22eY2kd77GgnwaBYKjNgD9n7ED+ccFVaaqRMZ13w7hBrR5fMmeMdWI8qM2gLyKrROSkiLwXU/5fRGSviOwWke9FlT8iIvudY7dElU8XkV3OsadFRFL7VrxrxaYDCY91L8xjyZwxPP/lq9NYIwOwaNEigM9b226/ig21/Nuuz2btCOFFAfv1KOBPJg5h7QPXsuzWCZmroEkomZz+c8A/A89HCkTkj4DbgM+p6nkRGeyUTwQWAJOAYcCrIjJOVZuAHwKLgSpgAzAX2Ji6t+Jd2w69H1eWL/CVL4yxD0YG3Xffffz4xz/eF11mbbttbmvsfL6sL+u/OjtDNTLt0WZPX1XfAGKj1gNAhaqed86J7HZ8G7BaVc+r6iFgPzBTREqBPqq6VVWV8B+Q+al6E152z7PVfHDuQlz5miXWE8q0OXPmAFyMKba23Qa3NXbeO/pRBmpiOqKjOf1xwBdEpFpENonIDKd8OHAk6rx6p2y48zi23JWILBaRGhGpOXWq7cEir6qsruONfafjyp+w7eG8rMvadja068rqOupDZ+PKJw/rk4HamI7oaNAvAEqAWcB/A9Y4eUy3XKa2Uu5KVVeqarmqlg8a1PpgkZf94LV9cWUThhbb4Ja3dVnb9nu7jiy3ELv94dhBvSy14yMdDfr1wDoN2wY0AwOd8hFR55UBR53yMpfyrBUIhqj/4NO48u/ePiUDtTHtYG07gUTLLcy3Hd18paNBfz1wI4CIjAO6AaeBl4EFIlIkIqOBscA2VT0GNIjILKfXdA/wUqdr72HfWLMzrmx4SQ9L63ifte0EEs25nzVmQJprYjqjzdk7IvJT4AZgoIjUA98GVgGrnKlujcC9ziDWbhFZA+whPED2oDO7AcIDZM8BPQjPbMjq2Q2Hz8TnPR+84YoM1MYksnDhQoArAbG23brK6joe/7c9LcqGFhfx/b+cbh0ZnxHVVlb/8oDy8nKtqanJdDXa5crHNvLphZaJzwlDi9m4dE6GamQSEZGAqqZ94wI/tetESyc/cfsUG5/yqNbatd2Rm2L3PFsdF/DBcvnGvxLl8pPd8tN4iwX9FKrYUOs6RXP+1GH2Fdj4UiAYos4lVQm2ro5f2SqbKbQmUB9XZmvqGL8KBEPctXwLTTEZ4B6FeTz2p5MsteNT1tNPkUAwxPufNLYoG1pcZGvqGN+qOngmLuAD3DGtzAK+j1nQTxG3RdU+N6JfBmpiTGokmorp7akfpi0W9FMgEAzxmz0n4soHFhdloDbGpMa/bj0cV5YncKfdjOVrltNPgSc31rqW24fD+E0gGGLt9np2BEPUHm+5tWdZv+7808JpNinB5yzop8C79R/GlS2ZM8Y+HMZXAsEQd6/YEre2ToQF/Oxg6Z1OqthQy6cxn5K+PQps2WTjO19/cWfCgA+wN6bnb/zJgn4nBIKhuM0kAB6eawHf+MvS1TsIvu8+Hz/CbsbKDhb0OygQDPGXz1TFldvSycZvKqvreKmNPZzBbsbKFpbTb6fK6jqeemUvpz9udD1uyy0YP0m0rk60fj0LeeiWK60zkyUs6LfDPc9Wuy6zEGGDt8ZPKqvreOyl99o8b8rwvhbws4gF/SS1FfD79yy0wVvjG8n08CMsrZNdLOgnYenqHa0GfIBn7p3R6nFjvKS1Qdn5U4fxSWMTJz/6lLtnXGa9/CxjQb8NFRtqWd/GIJdtdG78Zt7kUt6M6cgUFeTxV9eOsm+sWc5m77QiEAyxwmVKZjTbSML4TSAYcl1iwQJ+brCg34oVmw60uriUBXzjN5HlkmOXWAD41e7jGaiRSTcL+gkEgiF+Wxu/iFrEjFElFvCN7yRaLhlg7qSh6a2MyQgL+glUHTxDc4IPhwDL5tnXYOM/JT27uZbPnzrMUjs5os2gLyKrROSkiMRN6BWRvxURFZGBUWWPiMh+EdkrIrdElU8XkV3OsadFRFL3NlKvpGe3hKmdx23gNissWrQI4PO51LZDZ91vKhw7pDjNNTGZkkxP/zlgbmyhiIwAbgbqosomAguASc5zfiAi+c7hHwKLgbHOT9xresnre0+6ll8xqJeldbLEfffdB7Avtjxb23ZldR3//Nu4t0t+niTcMMVknzaDvqq+Abzvcuj/AA/RciOd24DVqnpeVQ8B+4GZIlIK9FHVraqqwPPA/E7XvosEgiFeTZDPXzR7TJprY7rKnDlzAC66HMq6th25GetczDKaowb0ZM3919g31xzSoZy+iHwR+IOqvhNzaDhwJOr3eqdsuPM4tjzR6y8WkRoRqTl16lRHqtgpKzYdcM3n2+Bt9uvKtp3Jdp3oZqwR/XtawM8x7Q76ItIT+CbwLbfDLmXaSrkrVV2pquWqWj5o0KD2VrFTAsEQr7hsfQg2eJvturptZ6pdV2yopergGddjtsRC7unIHbmXA6OBd5zxqjJgu4jMJNzLGRF1bhlw1Ckvcyn3nMfW73L9xM4YVWI9ouyXdW27YkOt654P3fKF73xxsn1zzUHt7umr6i5VHayqo1R1FOFGP01VjwMvAwtEpEhERhMe1NqmqseABhGZ5cxsuAd4KXVvIzUqNtSy55j77kDWy89+2di2E91wdfWYARbwc1QyUzZ/CmwFxotIvYh8OdG5qrobWAPsAX4FPKiqTc7hB4AfER4AOwBs7GTdUyrRLlhgSyZnq4ULFwJcSZa27crqOk59fN71mKV1cpeEJxx4V3l5udbU1HT5dRY/X8NvXHL5S+aMsZtWspiIBFS1PN3X7ep2nWjp5H49Cnlorm2Iku1aa9d2Ry6JB28nDC22gG98KdFsnZmj+1vAz3EW9IEnN9a6Dt7a1ofGjwLBEKcb3NM6Jz76NM21MV6T8+vpV2yoZdvhUFz5zROHWB7f+E4gGOKuFVtoanY/fvcM6+XnupwO+oFgiBVvJhi8vf7yNNfGmM5bt73eNeAXFeTx7f84yVI7JrfTO+u21+M2jj1n7EDr5RtfSjQtY0RJDwv4BsjxoL/vRPyc/JH9e/L8l6/OQG2M6bw7p5WR73KPsK0ZZSJyNr0TCIZcc/nXjR3ocrYx3hcIhli+6QAlvbtx7nwTBfnCoOLuLLputPXyzSU5G/RXbDrgWn7ntDLXcmO8LBAMcfeKLcQsosnDcy3gm5ZyMugnmpdv6+sYv6nYUMuamiOcbWyKC/gAqzYfsqBvWvD8HbkicgoIpuS1uvXohcjQ/J79yCvq2a/FQYULoT/8ThvPfZKKa3XSQOB0pivRhmyp40hVTe9SrqSmXef3GTQ8v2ffSxvbNp39kPyefVucoxcvfHrhdHB3Z66TYl5vN16vHyRXx7HAVlWN3wDL60E/1USkJhO33beH1TE1/FDHVPLD+/V6Hb1eP+h8HXN69o4xxuQaC/rGGJNDcjHor8x0BZJgdUwNP9Qxlfzwfr1eR6/XDzpZx5zL6RtjTC7LxZ6+McbkLAv6xhiTQ7Iu6IvIKhE5KSLvRZX1F5FXRGSf829J1LFHRGS/iOwVkVvSUL8RIvKaiNSKyG4R+ZoH69hdRLaJyDtOHf/ea3V0rpkvIjtE5JderF8qeb1dO9f0dNv2S7t2rtt1bVtVs+oHmANMA96LKvsesMx5vAx40nk8EXgHKAJGE97fNL+L61dKeLNtgGLg9049vFRHAXo7jwuBamCWl+roXPfrQCXwS6/9f861du2Htu2Xdt3VbTvjjbmL/oONivlw7AVKoxrmXufxI8AjUef9GrgmzXV9CbjZq3UEegLbgau9VEegDPgtcGPUB8Mz9eui9+ybdu1c17Nt26vt2rlOl7btrEvvJDBEVY8BOP8OdsqHA0eizqt3ytJCREYBVxHucXiqjs7Xy53ASeAVVfVaHf8ReAiIXnHGS/VLB8++X6+2bR+0a+jitp0rQT8Rl5XHE+5DkdoLi/QG1gJLVfWj1k51KevyOqpqk6pOJdzrmCkik1s5Pa11FJE/BU6qaiDZp7iUZfNc5Yy+Xy+3bS+3a0hP286VoH9CREoBnH9POuX1wIio88qAo11dGREpJPyheEFV13mxjhGq+gHwOjDXQ3W8DviiiBwGVgM3ishPPFS/dPHc+/VL2/Zou4Y0tO1cCfovA/c6j+8lnGuMlC8QkSIRGU14ZbptXVkRERHgWaBWVZ/yaB0HiUg/53EP4Cbgd16po6o+oqplqjoKWAD8u6r+pVfql0aeer9eb9teb9eQpradjkGTdP4APwWOARcI/xX8MjCA8MDIPuff/lHnf5PwiPdeYF4a6jeb8Nevd4Gdzs+tHqvj54AdTh3fA77llHumjlHXvYHPBrs8V79cadfONT3dtv3Urp1rd0nbtmUYjDEmh+RKescYYwwW9I0xJqdY0DfGmBxiQd8YY3KIBX1jjMkhFvSNMSaHWNA3xpgc8v8BE02iXRBuwx4AAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 432x288 with 2 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "fit,ax = plt.subplots(1,2)\n",
    "ax[0].scatter(ATL06_gdf_sub.index, \n",
    "              ATL06_gdf_sub['h_li'],\n",
    "             marker='.')\n",
    "\n",
    "ax[1].scatter(ATL06_gdf_sub.index, \n",
    "              ATL06_gdf_sub['srtm_height'],\n",
    "             marker='.')\n",
    "\n",
    "ax[0].set_ylim(1300,2800)\n",
    "ax[1].set_ylim(1300,2800)\n",
    "\n",
    "ax[0].set_title('h_li')\n",
    "ax[1].set_title('srtm_height')\n",
    ";"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 49,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXwAAAD4CAYAAADvsV2wAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjAsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+17YcXAAAgAElEQVR4nO2df5AdxXXvv+feXckoCLEWP/RjpRUyWHG05Id2DetyEkNwfpBgCMIYjPOe856NnCqSF1eSisEkKp6qcDkv5RSpelQSmeSV84Jk0A8sWxXy+GFIjMsrtKsYS0IWyAuSFv0WKyEiot29t98fd/pu397umZ6Znjtz955PlUq7d+fOdE93nz59zunTJIQAwzAMM/Mp5V0AhmEYpjmwwGcYhmkTWOAzDMO0CSzwGYZh2gQW+AzDMG1CR94FULnsssvEsmXL8i4GwzBMSzE8PHxKCHF51HWFEvjLli3D0NBQ3sVgGIZpKYjooMt1bNJhGIZpE1jgMwzDtAks8BmGYdoEFvgMwzBtAgt8hmGYNoEFPsMwTJvAAp9hmMIwfHAMj75wAMMHx/IuyoykUHH4DMO0L8MHx/CZxwYxPlnFrI4SHv/8APp6uvIu1oyCNXyGYQrB4MhpjE9WURXAxGQVgyOn8y7SjIMFPsMwhWBg+XzM6iihTEBnRwkDy+fnXaQZB5t0GIYpBH09XXj88wMYHDmNgeXz2ZyTASzwGYYpDH09XSzoM8SLSYeILiWizUT0YyLaR0QfIaL3E9GzRPR68D+3IsMwTI74suH/NYB/EUL8NICfA7APwP0AnhdCXAPg+eB3hmFSwqGLTFJSm3SI6BIAvwzgdwFACDEOYJyIbgNwQ3DZNwC8COBLaZ/HMO0Mhy4yafCh4S8HcBLA/yGifyeix4jopwBcKYQ4CgDB/1d4eBbDtDUcusikwYfA7wCwCsDfCCF+AcB/IIb5hojWENEQEQ2dPHnSQ3EYZubCoYtMGkgIke4GRAsADAohlgW//xJqAv9qADcIIY4S0UIALwohVoTdq7+/X/CJVwwTzvDBMQ5dZBogomEhRH/Udalt+EKIY0R0mIhWCCH2A7gJwKvBv88C+Grw/7a0z2IYhkMXmeT4isP/AwCPE9EsACMA/htq5qIniehzAA4BuNPTsxiGYZgEeBH4QogfAjAtJ27ycX+GYRgmPZxLh2EYpk1ggc8wMeBNT0wrw7l0GMYR3vTEtDqs4TOMI7zpiWl1WOAzjCO86YlpddikwzCO9PV0Ye0tK/H0nqO4uXchm3OYloMFPsM4MnxwDOu278X4ZBU733wbKxbMZaHPtBRs0mEYR9iGzySlKNFdrOEzjCPShj8xWWUbPuNMkaK7WOAzjCN85iqTBNPKkAU+wxQMU1ZKTlzGxKVIK0MW+AxjoEjLcKa1KdLKkAU+wxgo0jKcaX2KsjLkKB2GMcCbrJiZCGv4DGOgSMtwhvEFC3yGsVCEZTgfZ8j4hAU+w4SQp8BlxzHjGxb4DGMhb4HLjmPGN+y0ZRgLeadSYMcx4xvW8BnGQt4bZthxzPiGhBB5l6FOf3+/GBoayrsYzAwkqS3e9D12pDJFg4iGhRD9Udexhs/MeEy2eABOQluP1Mnbrs8waWCBz8x4dFv8ll2j2LprNJHQZkcq08qw05aZ8ejOTwKcnbF6HnN2pDKtDGv4zIxHd34CwJZdo5HOWJv5pgiOVPYjMElggc+0Bbot3kVo28w3ee/A3bDjENZu24OqEOxHYGLBAp9pS6SAlOYck8C0hWU2W7tWnwcAa7ftwWS1Fl03zn4EJgYs8Jm2xCXaxmS+aXaUjv681au6UalOhVKXiNiPwDjDAp9pS1yjbXTzTbOjdPTnEYDZnSWMT1RRKhHW3dbL2j3jDAt8pi2Js4tWNak0e/et/rzVq7qxelU3O2yZRPBOW6ZtcbHFp9m01exy8iTQvvBOW4aJwCXaxmTCue/Gq5sqVKPKybt/GVd44xXDBOibrIDW2GiVd1ZPpnVgDZ/xTiuaF4YPjuHT63+AiYpAZ5mwcc1H6pp1ETZahZF3Vk+mdWCBz3ilVc0LW3aNYrwSxLZXBLbsGq2XO++NVlG0wqTkm1ZUKoqAN4FPRGUAQwDeEkLcQkTvB/AEgGUA3gTwKSHEmP0OrQt3vilaLbmYbLtT5y40fE45lUfHtW8VfVLySasqFUXAp4b/hwD2Abgk+P1+AM8LIb5KRPcHv3/J4/MKAXe+RlrJvKC2XUeJ0FkmTFZEPfwxb7hvmWk1paJIeBH4RNQN4LcAPAzgj4KPbwNwQ/DzNwC8iBko8LnzNdJK5gW17SpVgbuuW4rFl140rdx5reDS9K1WWHUmLWMrKRUqRWgTXxr+IwD+FMBc5bMrhRBHAUAIcZSIrvD0rELRqp0vS4pmXrANNL3t7ljVPa3ceWrZavnKJcKRM+9h+OBY5PNbYWWQpoxZKxVZCOaitElqgU9EtwA4IYQYJqIbEnx/DYA1ALB06dK0xXHGV6O2kkbbjtgGmmz/tbesxNj5cWvb5bmCk31ry65RbB4excaXD2HLrtFIYdEKq860ZcxKqchKMBelTXxo+B8FcCsR/SaA9wG4hIj+CcBxIloYaPcLAZwwfVkIsR7AeqC209ZDeSLx3ahF02iZKWwx6q7t38wVnEkJ6evpwuDIaUxWptfBpmS0wqqzqGXMSjAXpb6pBb4Q4gEADwBAoOH/iRDid4joLwF8FsBXg/+3pX1WGtTBVJTZlvGDLiijct+Etb9+r2at4FQlpES1pGj3XF9b8ep16JozK3TCilPmvOzKRV0ZxxXMcaKoilDfLOPwvwrgSSL6HIBDAO7M8FmhyMF0YaJmC/38L16VeLZt5gDhHCrR6Ku1tbesxLrtexuEoWmg2fLcmwRpM1Zw6iRUFQJrt+3BigVzjZOOi8LiUua87cpFXBnHnSzjvL8i1NerwBdCvIhaNA6EEKcB3OTz/kkZHDmNCxNVCACTVYHHXnoD627rDbXdmnBtYB9C2OVZeQ/YvAhbrT2952hk7hvboI6z8vM90Q4sn48SEapBMsNqVTQ8XxcWPswDvNI14yqYW/H9tcVO24Hl81EuUf2UoKoQGDs/jvtuvDrWfVwa2JcQdnlWUTpcs1c9ukavCr+bexdi55tvRwpD06B2Xc6byhBXeTCVZ91tvbWjC6sCszrDy+7DPFAUu3JS8ljd5pkq2wdtIfAbBlNwDmiSxnFpYF9C2OVZRehwzV5l6O937Pz4NOG3YsHcaYLARTi4ClK1DOOTVW/ny95z/VJj2W1lTfueXeqbpVBNc+88VremZxbBLh+HGSnwTR0pzmCy4TJAfAlhl2cVwRHU7FWG6f3qwk//PY5wcBGkahmICJWqgICfzVHNtvOGPS9LoZr23nmsbouQKjstLSfwo7SCsI7kSyvy5Y2PqotLefN2BDV7lZFkktM18keeew1f/PgH63+LO1mqZeiaMwvrtu+NXX9bds4ikaVQTXvvPFa3RVhRp6WlBP6GHYewdtseVKoCszvNWoFPx1vSJWeWERKmMiUpZytvPIs7ycmBKvvFS6+fwo6R0wARJivJNEy1DElWjrbsnHHaJWsbdpYCLu298+p3ea+o09IyAn/44BjWbttTd7yOT5iFeVRHkoOka84sPPSdKc1s472NAz5rG2ESDcd23F7ccrbbxjM5UB957jW89PqpmvmlIgDUTDHjE9Faf1isf5L669k4CfHapRk27CwFnI9759Hvit7Xo2gZgT84croesgYApRIZtYKwjqQOEml7BWrL/K1K/nP5vCxthEk0HNuu0bjlLEp0jwu+VjR9PV344sc/WI/gKZeopuFPVlEF8P0DjVq/uvnJJdZfL0dUGVev6sam4dGGw8njtEuz2jBLASfvK/txUfvgTKJlBL66LJeD0bYEtnVSdZBAO7xdwC3kyjaQ4wqhJBqOrUxxJ45WsUVGrWj0HalR6O8cAB557jV8/8CpmuCsiEDnb9z85BLrH7U6BDBtRbDx3gFs3TUK2RPjtEurtGEY7bqPJE9aRuCbBGTcDqMOknK5hGq1ikoV6CwTehfNiwy5CkvElaTjxtWebJNE3ImjGbZIH/blqBWNKpTl9XF3JetafyVIlQxMbX4aWD4fHSXCREWgXKLIWP8tu0brG/0mJqvYsmsUW3eNGvvHluDzrUFStMc/3zgJ2Mrt0oZF34Wtt++WXaOFLu9MoGUEPjAlIOVh00fOvOecE0V+X9fwbDs2TSFXtmV0lstrk61Yv3eSZberYznJAPSludm0WH1Hqi5QTRuhwtImqH1i/7Fz5s1PRAAEQIQVC+aGmg03D08J7HKJQDCb3Uz9ZmD5/IZJIMx8FNaGraA9NyhgJcLm4dHETvR2ZvjgGMoXz1/gcm2hBb5J4KiROp1lQke5hErFLScKYB4k0ombdKOTy/I6aSRNXoM2zbN9TYA2LVbfkaoK1PEJ80aosDKpfaKvZ/rGrUdfOIDJSk1jr1TC468HR2qZLYGaI/bO/iVYvaobW3aNGvuNunIwKR6q+Uh1Luflp/G5aujr6cLaW1bi6T1HcVFnGc/tOx67vEVfxWSNHKfln+pa7HJ9YQW+zQ6qRupMVgTuvn7JtFOKbKYAk/MvzhZ5mwCymZvUlUSc6IuwVUeSQZBkUNie7XIvn/Zl0wStb6IDUBeoFGj/ernj+GT0Z6axra8ODlWRgu3m3oWN9VFWDqbvS/PR+MSUc3nnm2/HMl/6svH7VkCGD47VVy8dJbPy1szytCJynLoewlxYgW8T2nqkjumUIr2z29LJ6s9wya9jW0arn+sdcfWqbifBHZUnJspxbLtHVESJCZPACBtgeplcfARptDO9HaI2Qskybdk1Wh8brgIjjs/DNvnLNvjBT2r9+J7rl9ZXA9LWv3XXKB6+/dqGyUFObg3OZYujWDf9+fbTDI5MJSG0hUXHvZ8cF7UjJqcrb67fL3q0WVbIcQqhRaFYKKzAD4tI0SN1dPTObusYWUU66M8juEXSmCaguI7quBElNkwC49EXDli1fpt93EbcmHNXwRVmXwdQt/Vv2TXqPBHL9+EqTPRr1TZRHc0Dy+ejo1zrzwLApqHDWLloXn1y2Pnm2/UUyapzWe9DzUrr3DVnVt03UQ1+T4M+/kzKW1R5SsEKqVUjldIix+n1j5w54nJ9YQW+TUMxfeayLDcJ3KRaUJQAMi3rZZx1XFNImPAwCSmbSSDJpGYybXSUp6KcVMd3VnsBkqaJNq3Ukk7EaRlYbk59fN+NV+OTfd3YuONQzT9QFdYJOqyvNsteP3Z+HCUCqgIoETB2fjzV/dOsQuSqqVKt+T/W3rKy8Np9Vv6Gvp4uVN49fczl2kIJ/AMn3sWafxzCFz72AauGon/mkpMkrGPF1YJcBFDU4JTXyPvFMYVErUpM90ibNK4BuXJUVpBJVkqu39mqhTiaIrHCorVMzxyfqNn6Vy6a5zQRJ0FvV1vq4ztWdWOr4tANm6BtfbVZ9nqfz1HfT9w05cDUJFfbUiNSTz5ZUxR/Q6EE/nsTFTzz6nF8d/8JPBEI7qhZ0ZaTRCdMsOsO1rDnuWpTponJNS1CXBtx1LN9Le0HR05jMsgMWVEO6Eiiqbl8Z/jgGDYNHZ4KcSybI7FcHX59PTXnqYziWbd9r3U14Ct1r7o5zDTxyjLpNnvf7zJO2fVJVE045+M5Ps4TyMocmxVF8TcUSuBLJis1YSJjosOSpZlyksRBFxxRCbWSdjR9M47U9JPa1/PoLGF1T1KmqO/ICQaotesn+6ZsvEkdfmPnx+tRPOOWd+4zda++OezImfewZddovf6qQ1e12ft+ly6YJlF5ePpLr09FByXRyFXU92MLo43C5yTXDIoyQRVS4HeUCV1zZjklSzPlJIlC1d4aTAZKQi1bOGeSjmbajKM6oaM6QVFijZs9yExOvbC/6SsqUzm75sxC0KVQFbXf9WtNk3Ocupps9lt2jWLzcM1hDACbhw5j45qPNAi/CxPTczo1E9Mkevjt81MJ5zxppmrb2cJoXVAnf/X3IlKUCapQAv+izjJ+7WeuxBc+9gEMjtiTpekDdOO97i+yQYsJ0itMmQwIJaK6acAWzhlXm5Lhd8DUZhz5fRezhi2dQx6dp5mriyjfi0uSPF1rHDs/jiDyHSUAe46cxbrtexsOuLdNznHKbdocNhEIe6CmXMiy65E6q2NGq/hiYHnjRjA5wSZ1+ttQ2y7peQJAceziruS1MlcplMC/+oqLsf6/9td/tyVLi2v3VtFtaRIC8Kn+JbhDceLZ7G5xha0pakfiYtYw7UdIkhI5aTREnlpJ2Pux/S3MXjqwfD5md061BQENB9yv/95I/T765BwH0+awTYqG31mm+jvVI3VyjSfXNoKZJlYffUJtu6RBBUWxi7cShRL458crePCp3RCoRS+YNLi0jawK33K5BAhRS9OgmAXU+8XZfGQjzXLOZPuL+w7SHLZSpNWFK1G+Bj2f0hM7D9dNh0LI/DfC2UQIuJ1eJrNjnjh3AVfMnV3/mx6pM7B8PjbsONTgyG0GgyOnMRGsNOReBX1Fm0artr2jpJpvUezirUShBP7IyXfx+I5DAKZsnLqDKG0jmwa8baDKa9XdmUknnKSd2jZZxHkHSVcqvlYXaUmaelptN/3v6n30A+7jRo2YonJWLJhrfU/y8y1Bhsy+nsZInf3HzuHLT+0GAHzv9VMAkFjox0mzoW6sAoAnhw5P843Y+kRU+0RNFEnPOCiCXbyVKJTAVzubtHHahHCaRtYHfNg99h87hyd3HkYlcLzZ0h34wmUTWdx3YJokTQNw/7FzDVqlj9WFj/eRdIJRd9WGfc81FNImPNUQRhmVc9eHlxjfky1aqyFS58q5Dc99es/RRAJff3e/+5FleOylN6xRb3ose8UwBvU+YfNz6YT1G1lO6UcxnXHge3Xgk6KveFUKJfBVbUzaOE3EbeQ09ms9UsiU7iDOc+PkwQkTUnHegWmC0NMk/N2//gTPvHocQKNWmWZ1kVVOfJd7xY20iXqfYTmKOkoEIqpvRqsGexWkD4qI6tFAJoewXscrL3kfgLP1Z9/cu9BaJtcV2vhEFeu/N1KPUDJFvdU3pml+Bv09uaQt0QlbmQ+OnG7wo8gwVh9mpKyJW7YsJoeWTY+8/PKLcfv1S+s2fKlFyQMhwnJt2F5k2jS/pkghF2EbZ6OV+rwwLcjnqkYfgMff+c+G621aZdjqImmG0DDqQmhiSnBGoQtWIDrvS1zzlpoCoVIV+JUPXYkXfnyiHpVzx6pu9C6aV99Hsm77XtyxqtsaraW2xRc+9gHcsOKKUBu+S7/Wwx/lwS6A+YhQk59h/7FzkStOFwXAZEp99IUDGFheO12uXKK6YiVX0778dlkSp2xJZFFUv5Tp4lsyPfKcWWU8fPu19d+HD47h018fnBa77KoZDx8cwyPPvZa4s6gaT1iyNhM2W2ecPDhqGKpv56k+APcfO4dXRnfX/35z78JYSbn0a++IkSE0ytmp745VtT/bu5eCFajF24d9L67wNKVA+L2PfQC/F4QTy7pIhUGuMqTWLwMGRPBs2RbqSVf3XL+0QdDr7ynKli5/lv4IGf7okngQQN3EIlDLm5MkjYjp3rLf6u973W29+PNv7UZF1Ey7m4dH6wpekZ2zccrmO9iiwQLR6umRgamoAYnNrh/mXFQ7bdzNTa7abJwlbNw8OFH1S6MtqIJb/q9qlbbsmCb0MqrCLexQGJct9uruWJeVj3z3su2jzDquTm2XFAj6BK7HtcuInE1Dh/HNlw/VjzYEph93GGbSCLOl23aMu4Y/yvdRz4wZM41IFPr73rJrFIsvvQg3fehKPPvq8YaDZuS9w5zweWIbs2H9Mm2whfr3qltW5DqFE/iqCad30Tx0RtgUgfDQRbnB5qNXX2Y9KShK8CSxtds6Qtw8OFH1kzs0bTmE4pRZ1yrjdFD9Winc4tiYpflDd9zFXfmoGvOmocP1sFtb+V2c2qrNPlYKBENc++BILWWE6+rPNPDvu/Fquy3dsGPctjILex/y0BWTspRmdam+b/Vow44SobPDnhPJ1QmflOGDY/VJJc7mN/29RvXLNMEWpr+PT1ZbMx/++fFKgwlnVpnw0K292Hvk7DQbvkuWSfVlhR0LZxI8Ybk9XJdmpgFmm0CiJgFT/TpKhPGKmLYEjqqjq2kr7nLddG3Yd0w2Zum4+7Nv1cxL91y/dJp27WLble/ZJRumqexffmp3g9PX5VwBk9lFPxqxr6drKo+7EE6rP9vA1/uSKkSh7BiPawJR30fXnFnO5wMnuf+RM+9h48uHInMiZW3H37DjUN2kBNQ2yW28N9mk4tIvXYgaf+rf/+Cvjux3uWehBP67FyZR1kw4Y+fHG+z6QPgMKokjrHTBE5XbI61N0dW5aTLBqPW7s38JNgQ7NOXS2FbPpGW2dVCX8NGo7+iCRY2Iqgrgwad248X9J3DDiiuM2rVLnVwHmHqdKZom6lwB13TCMllaVQiUtDzuSdJI2K4BomPjTbikLU4rfPX+r573a1Nawpz3aQMapC28oujIaSYVnz6HqP4r//774+/9h8v9CiXwL57dgUkHE04aDduELnge+s7eumNNNR9ErShccHVuumhRq1d1Y9PQYSctP22Z1e8ljTZwcQD/2bd210MHBYBnXj2O54PoF5OJIouNN6rTV0bTRMXqh5ldVGeseh2hMY97WH916csmpSAOru2aRqCZnuHShnKVpzvvgfTRYIMjpxsimIBwf18UWfVLHxRK4M+ZVcbf3TsQGYaZtMOFaQJysAwfHJt2yEfa/D3q8x957rWGc0H3vHUWHSWaZmd2mdR0LV81G5hwKbOLcA87lMSGy3ekzV4V+kAtrr1cIggx3RafpB2i0PuXTK8Q9qywPrlp6DAmKgKbhw7joVt7U+1j8BHHbbuH3ue27hpNvOKwPdN0WM19N17tdA+T8x5AgzlW5u2P824GltfyK9VWD8BNH7qyfghTUuKsjJtJoQQ+4Gbjdu1wSeLCpUNNoGZSkhEEae2HquAUqGmOVQC73zqLjnIJd123pGGCs5kD9DqvXtVdi+6IEafuUkb5jkxCwHYoSdh9Xb8jhb5qT+3sKOGhT6zEniNnmxKlEUegqW1i+o5+QM+eI2ed+27cfRxhZXOxvzc4U8slbBqq5RgyPSvORKs+0/WwGhO2SVV1MH//wFTe/jhjdPWq7tjO2rgUYQNZ4QS+RHZWU/ywdORFdfYkceFdc2bVhYo0kzz0ifTpFFTBWSLgirmzceydCzVnVaWKxZdeFKpFAfbTseLGqYeVUV19yGerdZdOVWD6oSRh943zHWk+USMmgKn0A1lFaai4roZMdnupefb1dBkP6JH3Va/TMa3wAPcDc2zCJcqhKPvcW2fewzcDZ6przLiLuSvOYTU6ton48c8P4JHnXsP3D5yKrZTp70lNlhe2mTOJlp6149mFQgp8uXtMhujJ6A319KCo5ZL+cl3jwtdt39vgvKlU4qVTMJWlr6drmvZ06t0L9WttGq8qdMJi4sPi1OOgJs+qBr+bJp6tmpMtClPIZhS6wI2zJ0Anq2W03se27Bqthw5KIWs6oCfJRq+oSJ6ossn3FWUOVU2bsp3LJcKRM+/VN4npRNXH1P5J28E0Eff1dOGLH/9gorz9thQcYZsdk2rpPp25SUkt8IloCYB/BLAANTmxXgjx10T0fgBPAFgG4E0AnxJCjEXdT89fU6kKKGlKUFXyhYeZbPQkZ3HiwrX61QWfi+korFPo2hPQqPGqqxo9FC6ss4RFMESVVWXs/DhKhPoqxHQwdBL7rQ8n1sDy2kEhujM9iiyX0XqbqLn15QrpvhuvnnZAj8vkJVdueihqksgzPZTT1UkqNzttHh7FxpcPWVdWUZqrj/ZXibtJMuw+tsNubHVKo6VHlTGuYqLLPxd8aPiTAP5YCLGLiOYCGCaiZwH8LoDnhRBfJaL7AdwP4EtRNxscadw9Jk8heuylN+p5SqQ9O8xkY9PKXQcJESBAoWYSuVFDbhyx2b3VqBJde5KTkVofuVN6dmfjJjDbbkNXs46rJibNZzLhl82UNHxwrJ4PJaqDxrH5WlGc6aYcLzq6k9z3MlofwPuPnZu2QpLXhWm8YatN00avtJNsnHvIiKUwAVffVwD7Bjcv7Q+/E7gpGitKwUqrpdveQ9x66dfTrIt+yuX5qQW+EOIogKPBz+eIaB+AxQBuA3BDcNk3ALwIB4GvCx1ps//VlQtCtSSTySZuJ1MHibohxNTR5QuXwgSYuq7BfGNYDpsGo1ofAHUtUc/Tbttt6GLWUSei8cnpEQ3qxKEm/EoaNmoijf1TOtMnKyJyc5zuJC8hXaidDbWPDY6cjlwhye9EaaM+7L0+hGyUgKubQQPzq7qvIAts78W1P6r9zxaNBYRvJMwqFDhOe6vXX5ioonzRJU4d26sNn4iWAfgFADsAXBlMBhBCHCWiKyzfWQNgDQAsXbo09EVH2QWjTDYuqNqruiFE7+iq7Q+oaQjqRBO1HLbVp+7YRS2joSrE1d2eejoFF81Df8ZLr5/CjpHTuLN/ST06QU4cUiO2+T6SCKQ0SeAaV1/Rm+N0J3lYag1fxNH+ooSx671c3l0aU0GUgJPvWaC2u982yfnC9l5c+mPcPQC2NjJ9ntZPFHflIE2cdYXmfRdf5vIcbwKfiC4GsAXAF4UQ7xC5BdAJIdYDWA8A/f39AkgXmpnVcl3X7lXbX0cJuOvDSxtCulyXwyoyNGzlonkNGQ5lJ7i5dyF2jJw2brRSJxnbm5fXPPLca3jp9VO1VURFYMOOqQlpYLk54Zf+Hmwd1DVaQ408cdHM1PYIO/ha9YO4ptbwhU/tz+VettBNF3+S/L5r2KbLRNwMR6TtvbiUw9T/ovYAmPxzru8sziQQt+/09TSeheyKF4FPRJ2oCfvHhRBbg4+PE9HCQLtfCOBE0vurLzROaGZabPfXbX93fXjptPQPQDwtTQ8Nk8/Vd3fuPXI2dKNVVIKpvp6piAZTNsmB5fONCb9M99E7aNxoDVM0VdikqJbDtOt1+OAYPr3+B5ioCHQGeZjiHFXoA599MupeLlFCcU0gcVduroIqLBVo1bwAABpwSURBVCDBx3txKUfcySksgV7UOwPi7wCO23fUs5ABN7nvI0qHAPw9gH1CiL9S/vRtAJ8F8NXg/21Jn6G+0KjQTJ2kS62w75lsf2kiB8IGmd4JwjZauQ5W3eSkboKRk5k+objkzkkareErx49pk9NXDBNx1jRrN6UpSkh//zYhZ9v5nERjjxJUur8rLL++/r240WA+Jif5bP0sDVsCvbqpVBmTPnwwLvVVkqe95vIdHxr+RwH8FwC7ieiHwWdfRk3QP0lEnwNwCMCdcW6qO1dKgd0WaAzNDPsekCzPhmlm19Mlu2yKcu2wUU5elb6eRsfq2m17AKDhDNrxyejwTFkPk7nGtMPXV44VfVDGXcqGvdNT5y40/J5H7vQsw0CjzJoAjH6nO1Z1TztFzrbz2adZSqLa+YHo/Pqyrlm8Rxct2jRBmQ69Ud+ZHiWX1dnXpj7Q19PE5GlCiJdgH1s3JbmnqbHX3dZbe6FKaGbU91x31+o0RLMokTId5RI+2dfdYDcHzJuCAPfJRte4dSev3shj58cbUgmrKx49yiZqJeQigMPixuM4+FzLYCNsIgaAF/dPWQ07y9QQddEsstLsXJPPhSkicsObjHiSfGhB46HppvZI4xxWNWBbfn0d3+8xzmpBnaBKAK5dPA+9i+dhxYK5DX4kOc7lmFQDCZJs1nSpQ9pJsJA7bW3OlagTe/TvueyuNWGLCBmfrGLjjkPTTiTysaSTtkDdyQtMnzgGljeeAVoVUysePcomzCQTVpYo2zvgL6mcK6aJWIYD3vjTVzSkb1BjqptJVk7MOOa6MEVEN0HIfE6feWwwNJQxStCEXaMqAq42fJ/vMa6g1Ffc+46dw+63zjYENrikwvY9Fhr6f8JJsJAC39bYUS9Q/17SUE29g67bvjf0uDyTmeXzv3iV8+7XsHrbJr/6iieIRZfvSL9H15xZePCp3cZEWK6TgE1zb4ad0vZ+9ANTvvvjEw1ZR13SN2RBFiYRIJkADBtHcfLPuLSzi//G1V4P1ExTv3zN5bhs7uxUqRhcy6+itqG6H0eGQpuSKeonkGUxDrrmzKrv06kKJEqUWEiBn8Y04CtUU+2gKxbMjTwuTzezPPbSG1M7hLX83bZ66RON/N80aMPOVLVNVkBjNEdYNsSw9yFJqoUldWrKifXpPUexcuEleOylN+pavRACd354KRYlSMrlmyxWOUnGRNh3+nq6cHPvQvzgJ6cRtkMWiLfHI2r3sKksDdk0yyVUq1VMBhlOZqWcvIcPjuGVw2dAsG++CwtIkP4OGQr9xM7D+MTPLjTuLM5ydQvUZEwQP4cS7Bv7wiikwAfivzy10Wwn9bh8N2yZbDouTw03080se4++02DXM4XNAY0TgLyvy+HetnckP5dLenVzmDxDdGJy+g5hfRCq5fKZvyRseR3WFnq6AXVC1cNZZyJJBIrtO1E7ZJP4ZsJSDIe1u8kUK0mzcpRhujJyS9YTQD0lCBDua+vrqZ058fiOWu6rSlXgWz88AqB2xKh8b82IzBpYXsvbn8bMVViBH4c0zow439UHj/5dXQDpXn09bM42Aej5X8bOj8eexIDp+c0/2VcbkBtfPtQwCZg2L7nEH5veSRRhy+uottC/O/eiTjzxhY80JQRypqE6JvUdsnF9M6Z9JLbnmdpd76eqhp/Gfj84cjo40L1GtSqw98jZhr7sEtixelU3nth5uMHRDeW9+XCm2jZ0hUVlJenvM0Lgp7El+/zuOxcm8akPL2nQctS87isXzQuNm/67f/0JvvvjE3XTkLoETaJBmDqImjKirEUdqfVSsz66HODtStjSP6otkjjGmhUT32qEhQLHHRMu14e1u95PATSchZCmr5XLhMlA6AugtlM9ZmBHX09X3WdmGptpfVm2HdMuUVlxmRECP41H39d3palEZs1UtZzNSj70hz7RGEqoCt/n9x2v5+InAB+9ppb/BUh+bqfeQVy0BD0v/sqFl0yLP05jhzc9f/jgGN46817oaUhxNZwsY+JbHfkuTaHAcceE6fq42qmpn/pg5cJL8Mro2frvB07WwtVlaKhrYIfqMzNFGqWJKDJNGED0YTd5pUfOnTRLHV/ftWXXlCYboNaAe4+cbUjDIL//1pn3sDGwEwI1e6PM/5Lm8A9bucO+r+fFn3tRp9NGs6TPVw+86SwTbvrQlbh87uxEZVdpdhRRUYgTfWUKBY4bcWLS0LPQTsMw+ZxkGUxcu3ge1n5iym/hUq4wn1kaU4tpwtx/7Fxoyunc0iMXBdfO5JIiIMlzbdk19QQX+u/q97cGKRNKpVq+IJONsxkJqmymEwD1ycmHIB0+OIa//def4Ll9x+up7icqom7WemLn4XreJF/1mOmkiTlPE3GiXu9bQYnCVGd1slePFJX0Lp7ntUxpZYg+YUalnNZP6irNmjPXcOtpzBiB70KWS3zbLH/Hqm5sHjpcT+hlCzEL0xLU5bdruoA0tmub7V8NnVNj3pMI0g07DjUcVC4hgnUXcVx8OLlajTQx51HvKGmqkFcOn8GDT+3OLIrKVGd9IvsfN30QD31nKsPqykXznA/v8U2U0qlG15lSTuvZesslQnX8/DmXZ7eVwM96iW+a5ft6urBxjVskSZSWEJUJUxJnYrMNYr0s6rurVKq4+7rkMe/DB2vHWE4T9gBu/blF2P6jo8ZdxGH3c7UNz3SSrGpc3lHcaDapoDw5dBjPvHocALBpeBQb7/XvR7GtSHWtWUaprVw0zxp1ljUu7zGqDaUZDpjaVf6VZuXSaSXyWuLHFTomARZnsnK9Ns4g1t9dGm1tcOQ0KlXduFXjX/Yeq4e3VqqiftSiDXbM+slnFKW9q31qfHL6aWk60n81qczqWZp39ARxsgz66nRWR6kWfdYkk5P+Xl3GZlQbmsbiVxzL01YCvxWW+DYBFmeycr02ziTi893JDSS1tBPAVZdfjJGT7zbE16upI8KSwLWbYzbMOWmLmY8S5nG0TvmuX3r9FHa++XZo/p1NQ4cbPstCyTIlTNT/rqc5JpgzwvqeJE3v1XVshimJacZiWwl8wO8SP4sYb5sAi9PItmv18sZd8fh6d6bl9mceG5yWQyjqGEOgvRyzUc5J247pKGGua+9hWqd6WlpU/h11o9LPdTdGxbjWN6q/h9XfluZ49aruhl3zQPyos7jv1XfOnaRjse0EflJcNau0RG1QcX2GScszlTerFY9rugpJ0oNRTJNHXs64rHFxTprsvVErINekXH09U6elReXMOXLmvQbHfhJhn/YMBll3uVHKdK6xDKmOu0pMs9ksT78SC3wHkmhWSclKCIetHFwduq4rmiSToV6OuCsaWb5Pf31qpZCFgzBPXJyTcR2AQLykXGHPk6HJcgNiR7mEu65bkijbpev4CiuPXndV2DdEnZUodMOfCZf3WkQTMgt8B5JoVmnIQgOIW944+XR0fE2Gcd+Dvslt667RQgwyX9gESFp7r/SpxDHtAbV2lr/r5hOgFs21+NKLErVBnP5qq39Y3RuizqoCd123BID7KWmuwtxUtjzTfbDAdyCJZpUW350ibnl1oR0nn47pfTWjk0dtcpsJJFEGor4Tt2+ErXjDkvLFLXMWtm41u63aR3uVUM2osGfbvV1Im+gx7ftggR8QFcsdV7NKW5YkTqQk2obtu7rQtp3naXuOz1QMrrhucmOmY9LabUSteMslwp39S1JvtPI9vsKOyGxWtFfS5/jyGbLAh9vL9Nn54sQ8u3SKLNJDmya5qCMmTQmzgOZtte/rcd/kxjQic8fLyXLjmo84reBkpk3A7HgvEvq40tOONyPaK6kpeHCkMYNt0jHEAh/NjeW2pUJNEy6Zpvxh3zU5Um33lQnQ5FkAar1sp3aZ3k1e4WrtzpZdo/WDQsYrAltC/B9SGfjbIJ23nmnTZZXgSpo+ERWG3DVnVkNEVzMmrKTP0TPYJjneEGCBD6C5sdy6gDUdghK3U9jK7zJY0tRdtYeu3banHnc9bqiX7dQu9V7tvmM2T3RnpYvz8oUgyR0w5SSXSb3KQQLApInvgGxWrnJcnXtvYpqC0ixlIclz9Ay2SY43BFjgA2hu+JQuYPVDUFzCJV3K7zpYktZdvX8pOFBc5dS5C6HLZ5122zFbNFav6sYm5dwG06lVKoMjjekxSkQQQN3skDbxnXxGUjNGWBjy/mPn8LVn9tf3Htg2mxUJX0opC/yAZs7uptN90jakXv44AjRJ3dX7CyFAhHqKYyGAF/efiBXb3E47ZotIX08XNt4bb1Up02PIdN4rFszFk8pRgC6J78JIY8YIW/XWNPupa0tEhe9vvpRSFvgx8BVaqAvYLFYXWQtQef/xiSqqABAsNYWohUPK2OaT5y7gxDv/if3HznkNDWT8k3ZVCaAhB9KslP0ujRnDVr7pKxM0nD3RbOLIFB9KKQt8R7LOpe+7w8UVoHEnM3n/dd/Zi1dGz0IAIFHLzS1EbTv9JbM7sCE4xeuV0d0AEGrTZYdra2FqL3kUYJyzG2zEUVpcDzYyrUzS+BnSYJMpWe5ZYYHvSCvYmE1hkVnH+O47+k79546OxjN7H3nutYZrn95zNLfBxWSDTTi5nt0Qhqq0dM2ZZY3+SZKrvwgrSZNMAbLds8IC35E0kTDNII3QTjqZDY5MZUQk1A6YuOf6pfV3snLhJfje66fq169ceEmiujHFxNbnfCpH8nthfTvu84qykjTJlKwVSxb4jqSJhGkGaTpKUnu//r07VnVPeye//fOL8O1XjkAI4B++/wbOXZjM7Kg7prnY+pxv/9FW7fxWvW+3qsPfttrIsi4s8GPgEgkjP2+2xp+m0ydd5pq+p++qPf0fNUebQG1Dz4Ydh1It85niYOtzPs0m8iCV+vmt5el9u0hmGhNxjuDMui4s8FNg2rmXl8aftqMkXebq37Pl4JEaWtTBGUzrENbnfJlNBkdOY0I5JvGTfebVYVHMNDo+UoX7hAV+CvQOn7djtwid3iQEZNTG5uHRWDnHmeIT1ud8+LfUWHygtqFv+OBY7v3clbxlgg4L/JToHb4VbYm+MS1T+3q6cIdyrFyrDNiZQrODC3z5t9RYfAB49tXj+LfXT7aMSbBo/gUW+B4pui0RyDeqSN38ov7OZEsewQW+NFspMFvVJOgqE+RpYQRkGtTAAt8zRTCr2Mg7qijv57creZgVfGm2UmBu3TWKTUOHMVkRIKLE2SLzIEombNhxCH/+rd2QropNw6OZHc9Z8n5HDSL6DSLaT0QHiOj+rJ/HNDJ8cAyPvnCgrtmbooqahfr8CxO1TTlM9kjhW6Z0p1DFQQrqP/q1Fakn9r6eLjx8+7V46NZelIKD0ddu21Pfxd3KyNw+il+6vmEtCzLV8ImoDOBRAL8KYBTATiL6thDi1Syf286oJhsA0074ydOeOLB8PjpKhPGKgADwxM7D6F00j3ffZkxepkbfq92x8+OoVIW3bJy+SWIu1XP7SDYPjyY6/D2KrE061wE4IIQYAQAi+iaA2wCwwM8A3WRyx6ruaSmK8/Qx9PV04c7+JXg80MwqBRy0M5UimxpdGVg+H+US1Xd3V6rhB7U0k6TmSjW3j/RRALUD4FXTmy/fW9YmncUADiu/jwaf1SGiNUQ0RERDJ0+ezLg4xUE1tZh+T4JushHAtKV8X08X7rvx6twGyepV3egoTaXVkil0GSaKvp4urLutF+Wg+wjUNOE0Y8YXSc2lcvX1x7++Ag/ffi3e1znd9CYnk689sx+feWwwVX2z1vBNCfMa1i9CiPUA1gNAf3//9LXNDMR0mPK67XtTOzNNqQ6KFgopB62vFLpMe3HP9Uux98hZbNhxqJaGu1KMiJ20O91l+U3nRvt0umct8EcBLFF+7wZwJONnFh69AZ/ec9RLg9pstXkPBh2ZQrdIExHTOqxe1e3l0CCf+PKTmExvPmP5sxb4OwFcQ0RXAXgLwN0A7sn4mYXHln7AR4O2iq22VcrJFI8w4Zr3PhPXZ27YcQhP7zmKm3sXRgYt+HS6ZyrwhRCTRPT7AP4fgDKAfxBC7M3yma2ALf0Aa7xMq9MsgWsSrq2yz2PDjkP48lO1A4Fk+nAXoe+jLplvvBJC/DOAf876Oa2GLf3ATKQoZwYw2ZK3wC1a3hobT+85Ou13m8D3PXZ4py2TKXkLAaZ55C1w1XOWAeCVw2cKmWjt5t6FDQcD3dy70HhdFmMn8522THuT9+5epnnksaNXpa+nC2tvWQkioCKAZ149jk9/PV0YYxbcc/1SfOX2a/FL11yGr9x+rVW7z2LssIafAWzCmKJo2QKZ7ChC8sCx8+NQN64W1bRzz/VLI+32WYwdFvieKaIJI+/IhbyFANM88vZFDSyfj3KZMBkkp+koU8sqGVmMHRb4nsnbjqlThAkobyHAtBliSsVv9Z2cvscO2/A9k7cdU4dt6Ew7UUtGNvV7pcKpO1RYw/dM0UwYJjsg+xiYmcrA8vnoLNcysgLFULqKBAlRnEVPf3+/GBoayrsYM46wlMlF8DEwjE+adXpUkSCiYSFEf9R1rOG3Aaod8NEXDhTKx8AwvplpPiOfK3IW+G0Gh0kyTOvgO+iCBX6bUTQfA8MwdnxH/bHAb0Nm2pKXYWYqvlfkLPAZhmEKimlFnsamzwKfYRimwKgr8rQ2fd54xTAzDB/nIzPJyPrdp91IyRo+w8wgipBKo11pxrtPa9Nngc8wM4ii5XJqJ+K++yS2+LRRdizwGWYGwfss8iPOu0+zGkgTZccCn2FmELzPIj/ivPu8VmIs8BlmhsH7LPLD9d3ntRJjgc8wDNNk8lqJscBnGIbJgTxWYhyHzzAM0yawwGcYhmkTWOAzDMO0CSzwGYZh2gQW+AzDMG0CC3yGYZg2gQU+wzBMm8ACn2EYxiNFTk/NG68YhmE8UfT01KzhMwzDeCLuASXNXg2whs8wDOOJZqVITgoLfIZhGE8UPUUyC3yGYRiPxE2RPD5RBRGha86szMuWyoZPRH9JRD8moh8R0VNEdKnytweI6AAR7SeiX09fVIZhmJlDX08X1t6yEqUSoSoE1m3fiw07DmVq00+r4T8L4AEhxCQR/QWABwB8iYh+BsDdAFYCWATgOSL6oBCikvJ5DMMwM4ax8+OoCoGqAMYnqli7bQ+qQmRm00+l4QshnhFCTAa/DgLoDn6+DcA3hRAXhBBvADgA4Lo0z2IYhplpSLNOmVDX9F0jfJLg04b/3wE8Efy8GLUJQDIafDYNIloDYA0ALF261GNxGIZhio3q5O2aMwvrtu/N9NjDSIFPRM8BWGD404NCiG3BNQ8CmATwuPya4Xphur8QYj2A9QDQ399vvIZhGGamojp5VyyYm+mxh5ECXwjx8bC/E9FnAdwC4CYhhBTYowCWKJd1AziStJAMwzDtQNbHHqaN0vkNAF8CcKsQ4rzyp28DuJuIZhPRVQCuAfBymmcxDMMw6Uhrw//fAGYDeJaIAGBQCPF7Qoi9RPQkgFdRM/XcxxE6DMMw+ZJK4Ashrg7528MAHk5zf4ZhGMYfnDyNYRimTWCBzzAM0yawwGcYhmkTWOAzDMO0CSzwGYZh2gQW+AzDMG0CC3yGYZgC4/MYRD4AhWEYpqD4PgaRNXyGYZiCEvdQ9ChY4DMMwxQUNV++j5TJbNJhGIYpKHEORXeBBT7DMEyB8ZkymU06DMMwbQILfIZhmDaBBT7DMEybwAKfYRimTWCBzzAM0yawwGcYhmkTSAiRdxnqENFJAAdzLsZlAE7lXIa0zIQ6ADOjHlyH4jAT6mGrQ48Q4vKoLxdK4BcBIhoSQvTnXY40zIQ6ADOjHlyH4jAT6pG2DmzSYRiGaRNY4DMMw7QJLPCnsz7vAnhgJtQBmBn14DoUh5lQj1R1YBs+wzBMm8AaPsMwTJvAAp9hGKZNaCuBT0T/QEQniGiP8tn7iehZIno9+L9L+dsDRHSAiPYT0a/nU+rpWOrxEBG9RUQ/DP79pvK3wtWDiJYQ0QtEtI+I9hLRHwaft0x7hNSh1drifUT0MhG9EtTjfwaft1Jb2OrQUm0BAERUJqJ/J6Ltwe/+2kEI0Tb/APwygFUA9iif/S8A9wc/3w/gL4KffwbAKwBmA7gKwE8AlPOuQ0g9HgLwJ4ZrC1kPAAsBrAp+ngvgtaCsLdMeIXVotbYgABcHP3cC2AFgoMXawlaHlmqLoGx/BGADgO3B797aoa00fCHEvwF4W/v4NgDfCH7+BoDfVj7/phDighDiDQAHAFzXlIJGYKmHjULWQwhxVAixK/j5HIB9ABajhdojpA42ClcHABA13g1+7Qz+CbRWW9jqYKNwdQAAIuoG8FsAHlM+9tYObSXwLVwphDgK1AYwgCuCzxcDOKxcN4rwwVwEfp+IfhSYfOSyr/D1IKJlAH4BNa2sJdtDqwPQYm0RmBF+COAEgGeFEC3XFpY6AK3VFo8A+FMAVeUzb+3AAt8OGT4rcgzr3wD4AICfB3AUwNeCzwtdDyK6GMAWAF8UQrwTdqnhs0LUw1CHlmsLIURFCPHzALoBXEdEvSGXF7Ieljq0TFsQ0S0ATgghhl2/YvgstA4s8IHjRLQQAIL/TwSfjwJYolzXDeBIk8vmjBDieNDhqwC+jqmlXWHrQUSdqAnKx4UQW4OPW6o9THVoxbaQCCHOAHgRwG+gxdpCotahxdriowBuJaI3AXwTwK8Q0T/BYzuwwAe+DeCzwc+fBbBN+fxuIppNRFcBuAbAyzmUzwnZIQJuByAjeApZDyIiAH8PYJ8Q4q+UP7VMe9jq0IJtcTkRXRr8fBGAjwP4MVqrLYx1aKW2EEI8IIToFkIsA3A3gO8KIX4HPtshb490M/8B2Ijasm4CtdnxcwDmA3gewOvB/+9Xrn8QNc/3fgA3513+iHr8XwC7Afwo6AgLi1wPAL+I2vLzRwB+GPz7zVZqj5A6tFpb/CyAfw/KuwfA2uDzVmoLWx1aqi2Ust2AqSgdb+3AqRUYhmHaBDbpMAzDtAks8BmGYdoEFvgMwzBtAgt8hmGYNoEFPsMwTJvAAp9hGKZNYIHPMAzTJvx/1azNnQQLu0IAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "ATL06_gdf_sub['diff'].plot(marker='.', linestyle='none');"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 50,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "21.185592375756876"
      ]
     },
     "execution_count": 50,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "ATL06_gdf_sub['diff'].mean()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "#### geoid offset !"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.6.7"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}
